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I Abstract 

D . 

. We study the fast rotation limit for a Bose-Einstein condensate in a quadratic plus quartic con- 

fining potential within the frameworlt of the two dimensional Gross-Pitaevskii energy functional. As 
the rotation speed tends to infinity with a proper scaling of the other parameters in the model, a 
linear limit problem appears for which we are able to derive precise energy estimates. We prove that 
(~| ' the energy and density asymptotics of the problem can be obtained by minimizing a simplified one 

^ |. dimensional energy functional. In the case of a fixed coupling constant we also prove that a giant 

vortex state appears. It is an annulus with pure irrotational flow encircling a central low-density hole 
around which there is a macroscopic phase circulation. 



oo 
O 



1 Introduction 



1.1 Physical background 

, Since the first experimental achievement of a Bose-Einstein condensate by the Jila and MIT groups in 1995 



(2001 Nobel prize in physics attributed to Cornell, Wieman and Ketterle), these systems have been the 
subject of many studies from the condensed matter community. A reason (among others) for this is the 
fact that a Bose-Einstein condensate is a good system to study superfluidity issues, such as the existence 
of quantized vortices. These vortices can be generated by rotation of the container (generally a magnetic 
QQ ' trap) enclosing the condensate. They are the subject of an ever increasing number of experimental and 

, theoretical papers, see e.g. the review [35] or the monograph [T] for extensive lists of references. 

Recently, mathematical contributions have studied certain issues arising from the physics of rotating 
Bose-Einstein condensates. Let us cite some of these issues : 

X" 

^ ' • the appearance of the first vortices when increasing the rotational speed of the trap in the strong 

5^ , coupling regime [571 12H] ; 

• the formation of a vortex lattice characteristic of regimes when the rotational speed nearly deconfines 
the atoms [lll[5]; 

• the interaction of many condensates in an optical lattice [7]; 

• the symmetry breaking of the ground state of a condensate in parameter regimes |42] ; 

• the energy and density asymptotics of strongly interacting condensates in anharmonic traps [TSl 

dSKm . 



*Universite Paris 6, Laboratoire Jacques- Louis Lions, 175 rue du Chevaleret, 75013 Paris, France. email: 
rougerieOann. j ussieu.fr 
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Most of the available mathematical studies on Bose-Einstein condensates are made in the framework 
of the Gross-Pitaevskii energy. Here we will use a two-dimensional energy although the actual energy 
should be three-dimensional (see the discussion in Section 1.4 below). The consensate is described by a 
complex macroscopic wave-function ijj minimizing the energy given (in the rotating frame) by 

EgpW^ [ \\7^p-inx^ip\^dx+ [ {V{x)\i:\^ -n^\x\^\iP\^)dx + G [ \^p\'^dx, (1.1) 

under the mass constraint 

/ \iP\^dx = l. (1.2) 

JR2 

The matter density at some point x is given by |7/;p(x). A vortex is a zero of the wave-function ■0 around 
which there is a quantized phase circulation (ie topological degree or winding number). 
Here we note x — {xi,X2), x"*" = {—X2,xi), V{x) is the potential confining the atoms (generally repre- 
senting a magnetic trap), is the speed at which the trap is rotated around the axis perpendicular to 
the xi — X2 plane and G > is the coupling constant modelling the atomic interactions which we assume 
to be repulsive. The different terms in the energy represent the kinetic energy (including the effect of 
Coriolis forces due to the transformation to the rotating frame), the potential energy (including the effect 
of centrifugal forces) and the energy due to the interatomic interactions. 

The model based on the Gross-Pitaevskii energy is an approximation of the quantum mechanical many- 
body problem for TV bosons at zero temperature. We refer to [37] and references therein for a discussion 
of the connection between this type of energies and the iV-body problem. We remark that the Gross- 
Pitaevskii energy in the rotating case {Q ^ 0) was rigorously derived from the TV-body problem in |36j . 
but only for fixed rotation (ie fixed fi) and fixed interactions (ie fixed G). In the regimes we are going to 
consider (f2 and/or G going to infinity), the derivation of the Gross-Pitaevskii energy from the quantum 
mechanical many-body problem is still a mostly open problem. To the best of the author's knowledge, 
the only rigourous results available are those of 14j, where it is proved that the Gross-Pitaevskii energy 
is a correct first-order description in regimes similar to the ones we are going to study in this paper. 
We remark that our results provide an a posteriori criterium for the justification of the reduction to 
the Gross-Pitaevskii energy in some regime of parameters. However, it is far from a rigorous proof, see 
Section 1.4 below. 

In most experimental situations, the potential V is harmonic, of the form 

V{x) — aixf + a2X2- (1-3) 

Such a potential does not allow to take arbitrarily large rotation speeds, as the so-called effective potential 

V{x)-n''\x\^ (1.4) 

is not bounded below for fl"^ > min(ai,a2) and consequently neither is the energy (jl.ip . Physically, 
this corresponds to the fact that the centrifugal force overcomes the magnetic trapping force when fi^ > 
min(ai,a2), thus the condensate is no longer confined and the atoms fly apart. Fetter suggested [2T] to 
use instead a potential with a growth steeper than harmonic, for example of the form 

V{x) = \x\^ + k\x\'^. (1.5) 

The nice feature of this potential is that the centrifugal force is always compensated by the trap- 
ping force, and thus one can in theory take arbitrarily large rotation speeds. Experiments with this 
type of potentials have been realized by the ENS group [HI |33jj using a blue-detuned laser beam 
to create the quartic part of the potential. The experiments motivated numerous theoretical studies 
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[SlinillllMlllSllHISniinilSllSlISS] revealing the very rich vortex structure one can expect to be 
displayed by such systems. When the rotational speed of a condensate trapped by a potential of this kind 
is increased from zero, many phase transitions are expected to happen. Firstly, vortices are expected to 
appear one by one as is the case for a harmonically trapped condensate, but eventually with multiply 
quantized vortices becoming stable |29l 1301 [39] . which is never the case with a purely harmonic trap. 
When the centrifugal force begins to compensate the trapping force corresponding to the quadratic part 
in the potential, a triangular lattice of singly quantized vortices (Abrikosov lattice) similar to that which 
is observed in purely quadratic traps appears [Hj, but a new feature of the quadratic plus quartic trap 
is the existence of a critical speed for the centrifugal force to create a central hole in the condensate. 
The resulting state is an annular condensate with a vortex lattice encircling a central giant hole carrying 
a macroscopic phase circulation [6l [9l [191 US [32]. At even larger rotation speeds, a new transition is 
expected to happen, all the individual vortices present in the annulus retreating in the central hole and 
gathering in a single multiply quantized vortex at the center of the trap. This is what we will refer to 
as the giant vortex state: an annular condensate with pure irrotational flow encircling a central multi- 
quantized vortex [H [m 1311 [31 133] . 

In this paper we aim at justifying rigorously the appearance of the giant vortex state in the limit 
— !■ +00. The existing mathematical results on giant vortices in Bose-Einstein condensates (mainly 
those of [2ji to our knowledge) focused on the case when the trapping potential is taken so that the 
condensate has an annular shape even at slow rotation speeds and studied the effect of the central hole 
on the vortex structure in the annulus. Here we study the case where, as the rotational speed increases, 
a central hole is created in the condensate by the centrifugal force. 



1.2 Model 

Our model is the following: we consider a Gross-Pitaevskii energy of the form 

EGp{i^)^ f + k\x\^\M^ + G\yj\^)dx (1.6) 

to be minimized under the mass-constraint (II. 2p . In order to study the asymptotics of the problem when 
n +00, it is more convenient to change scales, setting for > 1 



We then have 
with 



u{x) = Ri) (Rx) , R = 



n'^ - 1 



2k 



2k 



VU — lll ; X U 

2k 



(I 



xl^-2\x\') \u\\ 



(1.7) 
(1.8) 
(1.9) 



Let us emphasize that formal calculations [23] suggest that the ground state of (|1.6p should be confined 
on an annulus of radius ~ i? so that our change of scales is quite natural in this setting. Because of the 
mass constraint 



^dx = 1, 



(1.10) 



we do not change the physics of the problem by adding to JF^ any multiple of Jjja Hence it is 

equivalent to minimize 



-^o(w) 



2k 



Vu — ifl- 



2k 



-X u 



' + G\u\'+^^^-^{\xf-l)'\uf 



(1.11) 
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Defining the parameters 



- 1 

n- — - 

2k 



Dn - (1.12) 



we have 

2k 

^n{u) - (1.13) 

with 

F^{u) = j ^ (^Vu - iujx^u\^ + Y (kP - 1)' |wp + Gl^l^j dx, (1.14) 

and we are going to study the asymptotics of F^^. Results for !Fq,{u) are straigthforward modifications of 
those we will present. 

The potential Dquj^ (I^P ~ l) is positive and has a degenerate minimum for |a;| — 1, so that we ex- 
pect the condensate to be tightly confined on an annulus centered on the circle |x| = 1 in the limit 
Dqoj^ +00. This is reminiscent of semi-classical studies of Hamiltonians with potential wells (see pS] 
and the references therein), although our analysis will be quite different. 



We denote 

Iu = mil^F^{u),ueH\R^)nL^{R^xUx), J \u\^ = 1^ = F^{u^). (1.15) 
It is classical that I^j is achieved by some (a priori non unique) u^j satisfying the Euler-Lagrange equation 

2 

- (V - iujx-^fu^ + {\^\^ - 1)^ + '2G\u^\'^u^ = n^u^, (1.16) 

where is the Lagrange multiplier associated with the mass constraint. We aim at studying the 
asymptotics of and Uuj when oj — > +oo. In particular we want to confirm rigorously the predicted 
transition of to a state with pure irrotational flow of the form 

(r,^)K^/(r)e™^ (1.17) 

where r, 9 arc the polar coordinates and n an integer. 

More precisely, it is to be expected that will be close to a function of the form (|1.17p with n ^ lo and 
/ a function tightly confined on a shrinking annulus centered on the potential well. 



1.3 Main results 

In order to state our main results, let us introduce some notation. For every n G Z we define the one 
dimensional energy Fn 

K(/) = 27r/ {\f'{r)\^ +Vn{r)\f{r)\^)rdr (1.18) 



defined for every / G H^{W^ ,rdr) H L'^iW^ ^rVn{r)dr). The potential is given by: 

2 2 2 

Vn{r) = '^-2nuj + {l^Dn)uo\^+DnY + DnYr\ (1.19) 
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The energies F„ appear when decomposing the function u^j as a Fourier series in the angular variable 
(see (|1.43|) below). We will denote gi^n the ground-state (it is unique up to a multiplicative constant of 
modulus 1, that we fix so as to have > 0) of the energy p.lSp satisfying the mass constraint 

2^ / 9ln{r)rdr = 1 (1.20) 

and 

Ai,n = -Pn(ffi,„) = inf |i^„(/), / e H\m.+ ,rdr) L\R+ ,rVn{r)dr),2TT £^ f{r)rdr = l| . (1.21) 
Similarly we introduce 

EM) = 27r / (|/'(r)|2 + Vn{r)\f{r)\' + G\f{r)\^) rdr (1.22) 

with 

7n = ^^n(*„) = inf |£;„(/), / e H\M.+ ,rdr)nL\R+,rVM)dr),2TT £^ f{r)rdr = l| . (1.23) 
Similarly to gi^n we can ensure the unicity of by requiring that > 0. 

We will denote and ^2 the first two cigenfunctions normalized in i^(R) of the harmonic oscillator 

+ (L24) 

associated to the eigenvalues 1 and 3. 

Our results concern the energy (|1.14p in the limit lo +oo. They hold when the interaction energy 
is small compared to the potential and kinetic terms so that the first order of the energy is a quadratic 
term, leading to a linear equation. We show that this condition is fulfilled if 

cj^+oo, t^>G2>.g2, Dne{Q,l), Dn^De{Q,l] (1.25) 

where D and g are fixed constants. The assumption that is larger than some constant is only a matter 
of convenience for writing the results. The opposite case where G <C 1 can be dealt with using the tools 
that we develop in this paper, and is actually easier than the case where G stays bounded below. We 
will refer to the limit p.25p as the extreme rotation regime. 

We have the following result (the notation a (x b has the usual meaning that a/b converges to some 
constant.): 

Theorem 1.1 (Energy and density asymptotics in the extreme rotation regime). 

Let Uoj be any solution of U.15\} and suppose that \1.25\) holds. There is an integer n* — w + 0(l) so that 



= F^{u^) = Ai,„. + 2^G / gi,„. (r)Vdr + O(G^) = 7„- + O(G^). (1.26) 

Moreover 

II - gin' IIl=(r^) < CGi/2 « \\gl^. |U.(R2) cx (1.27) 
and Iwojp converges to a Dirac delta function at \x\ — 1 in the weak sense of measures. 
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Note that |gi,„»P converges to a Dirac delta function at |a;| — 1 as a consequence of Theorem 12.11 
below, to which we refer for more details on Ai,„ and 51, In particular we have 

and 

f 4 {2Dn + Af^ f 4 
27r / gi,n'{r) j= uj ' / S,i{x)dx, 

thus Theorem 1 1 . 1 1 savs that the energy of the condensate is described to subleading order by a simplified 
one-dimensional problem, which also gives the density profile. The second equality in p.26p holds because 
the energies F„ and £'„ are actually very close in the parameter range p.25p , see Theorems 12.11 and 14.11 
below. Similarly (jl.27p could be stated with 5'„. replacing gi^n-- 

The energy En is obtained by restricting to wave functions of the form f{r)e™^, so Theorem 1 1.1 1 is a 
first step towards the understanding of the giant vortex. However, a result such as (|1.27p is not sufhcient 
to conclude that there are no vortices in the bulk of the condensate. Actually our method docs not give 
precise enough energy estimates to do so in the whole regime (jl.25p . 

A special case where we can prove much more detailed results is that of a fixed coupling constant: 

Theorem 1.2 (Refined asymptotics in the case of fixed G). 
Let Uuj be any solution of il.l5\) and suppose that there holds 

07^+00, G is a fixed constant, Dne{0,l), Do ^ £> G (0, 1]. (1-28) 

Then there is an integer n* — lu + 0{1) so that for any e > 

= FUu^) = 7„. + 0{uj-'/^+'). (1.29) 

Moreover there exists a G M such that, along some subsequence 

'*„-e™ «|U2(R2) <C,c^-i/i6+^ (1.30) 



Theorem ll.2l is still mainly concerned with energy and density asymptotics, although (|1.30p allows to 
identify a global limiting phase. We now state our result about the appearance of the giant vortex, that 
is an annular condensate with no vortices in the bulk : 

Theorem 1.3 (The giant vortex state in the case of fixed G). 
Under the assumptions of Theorem \1.2\ 



1. There is a constant a so that the following estimate holds true pointwise along some subsequence 
In particular 

- I*«*IIIl=o(e^) < ^0.3/16+^ « |lvl/„.|l^^(j,,) (1.32) 

2. |uaj(x)| uniformly on 

Tif, = |a; e M^ ||a;| - > Suj-'^ Iii{lu)^ (1.33) 

for every S > j^. 
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3. Suppose that there is some so that 

Uu^ixu.) (1.34) 
then, if u) is large enough, we have for every 5 < 

\\xu,\-l\^ >5uj~^\n{uj). (1.35) 

We remark that whereas there can be a great number (proportional to G) of integers n* satisfying 
p.26p and (|1.27p . there are at most two integers satisfying (|1.29p and p.30p . More precisely, two cases 
can occur : either 7„ is minimized at a single n and the whole sequence u^^ satisfies (11.301) and p.3ip 
with n* = n or 7„ is minimized at two integers h and n + 1. This case is very particular and happens 
only for special values of the parameters. For these particular values of the parameters, the sequence 
could in principle oscillate indefinitely between '^ne™^ and ^a+ie'^""'""'^^^, so that we have to extract a 
subsequence to state the result. 



We will discuss these results in further details in Section 1.4 and now present some ideas of the proofs. 
It is natural (see [22 for example) to construct a first trial function of the form 

vtest{r,e) = ft,st{r)e''' (1.36) 

V - r 



ftest{T) — CtestS, 



V 



where ^ : M ^ M is chosen with Jjj = 1, n = [w] is the closest integer to oj and rj a small parameter to 
be chosen later on. We want the mass constraint p.lOp to be satisfied, which implies 



cL, ~ (V27rry)-i (1.37) 



when r] under some conditions on ^, see Remark [T] in Section 2 below. 

Expanding the potential Dquj^ (I^^P ~ l) allows to compute the energy of Vtest- Minimizing it with 
respect to rj gives the main scales of the problem: taking 

r,= (2A2 + 4)-'/"^-i/2 (1.38) 

we get 

FUu..) < F^ivtest) < {2Dn + 4)'^' u J (C'(a;)' + x'^^^ix)) dx + 0{Glo^''') (1.39) 
so that when oj ^ the natural choice is to take for ^ the Gaussian f i achieving the infimum 



1= inf {C{x)^ +x^£_^{x))dx. (1.40) 

The symmetry of allows to improve the remainder term in (|1.39p : 

FM < F^vtest) < (2Dn + 4)i/'c^ + Gc.i/2(^^Zil±ilL / ^f{x)dx + 0{l). (1.41) 

V 27r jR 

This computation suggests that the dominant terms in the energy will be the kinetic and potential 
contributions, whose sum will be of the order of to. The interaction energy will appear only at second 
order, with a contribution proportional to Glo^/'^, which is much smaller than uj in the extreme rotation 
regime. 

Thus, it is natural to expand u^^ in a Fourier series in the angular variable 



.{r,e) = Y,f^{r)e^-'. (1.42) 
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Proving that converges to a function of the form (|1.17p then amounts to showing that /„ converges 
to for all n except one, and the crucial fact will be the following decoupling of the quadratic part of the 
energy (|l.lip : 



with the energies F„ being defined by (|1.18p . The proof of Theorem II. II is based on the detailed study of 
the ground states of these energies. In particular, for most values of the parameters, Ai_„ is minimized 
at a unique integer, n*. 

We prove that because of the properties of the potentials Vm some modes /„ carry too much energy to 
match our upper bound (|1.41|) and therefore can be eliminated. For the other ones, we use a lower bound 
on the gap between the ground-state energy of (|1.18p and the first excited level. This gap is of order ui, 
which is much larger than the interaction energy if w 3> , as shown by the preceding computation. We 
use this fact to prove that /„ oc when uj goes to infinity. 

This very particular form allows one to bound the quartic part of the energy from below. Combining 
the lower bound with the analysis of the one dimensional energies gives an estimate of the total energy 
contribution of each mode. Two effects come into play: concentration of the mass of u^j on the mode 
/„* cx gi.n* is favorable for the quadratic part of the energy because Ai^„» is the minimum of Ai.„ with 
respect to n, but it increases the quartic part. We then show, confirming the energy scaling of (jl.4ip 
that, in order to match the subleading order of the energy, it is more favorable for to have its mass 
concentrated on only one mode. 

However, the remainder terms that we obtain are not small enough to confirm this rigorously in the whole 
range of parameters (jl.25p . The energy expansion (jl.26p still allows for many modes (a number propor- 
tional to G) not to be asymptotically in the limit. We get the density asymptotics (|1.27p nevertheless 
because the moduli of these modes are very close to one another. 

If we fix G we can prove (|1.30p , which is much more precise than (jl.27p and requires a refinement of the 
method described above. In particular, using the Euler-Lagrange equation for u^, and the result p.27p 
we refine the estimate fr^ cx Qi^n: 

finding a second order correction. This in turn allows to refine the 
remainder terms in the energy expansions and prove p.30p . 

The dominant part of the Euler-Lagrange equation (terms corresponding to the variations of the kinetic 
and potential energies) is linear in the situation we consider. We are thus able, using elliptic estimates for 
the Ginzburg-Landau operator due to Lu and Pan |.38j , to get estimates in stronger norms. The uniform 
estimate (jl.Sip is obtained by this method. It implies that is confined to a shrinking annulus in which 
there are no zeroes, because so is 5'„. (see Theorem 14. ip . 

We remark that a Fourier expansion such as (|1.42p has been used in [8 to compute the lowest eigenvalue 
of the Ginzburg-Landau operator on a disc with an applied field going to infinity, but with no potential 
term in the energy. In that paper the minimizer gets confined close to the boundary of the disc, whereas 
in our work it is the potential that forces the confinement to a shrinking annulus. In [8], a rescaling 
similar to (jl.36p leads to an energy expansion in powers of the small characteristic length of the problem. 
We perform a similar analysis in Section 2, with the difference that our original problem is posed on the 
whole space M^, leading to rescaled eigenvalue problems on a line whereas in |^8J a compact domain is 
considered, so that after rescaling one gets eigenvalue problems on a half line. 



1.4 Discussion 

Theorems 11.21 and 11.31 confirm the expected qualitative features at high rotation speeds : the condensate 
is tightly confined on a shrinking annulus in which there are no vortices. We identify a simplified limiting 
profile with a Gaussian shape (see Theorems 12.11 and 14.11 for details on the functions gi^n* and re- 
spectively). We are also able to identify a global limiting phase and prove that all vortices gather in the 




(1.43) 
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central low-density hole, confirming that a giant vortex state appears when oj goes to infinity with G fixed. 



We prove that the width of the annulus where the mass is concentrated is of order w^^/^ ~ il^'^/^fc^/^ 
in the extreme rotation limit. This is larger than the width predicted in former physical studies 
and our limiting density is more regular than the Thomas-Fermi type profiles they suggest. This is due 
to the fact that these studies are made assuming that the interaction energy is the dominant part of the 
energy, leading to a profile solution of a nonlinear equation. Here we study a different regime of parame- 
ters and in particular we show that the Thomas- Fermi approximation for strong interactions breaks down 
in the extreme rotation limit, where u^j converges to the solution of some linear reduced problem. 

We used a two-dimensional Gross-Pitaevskii energy instead of the complete 3D energy. A rather 
natural question would then be to find in which situations the reduction to the two-dimensional model is 
justified. One can always assume that the condensate is tightly confined in the direction of the rotation 
axis, so that the problem is essentially 2D, but in such fast rotation regimes it is to be expected that the 
reduction is valid independently of the strength of the vertical confinement (see [1] where this is shown 
for the case of a fast rotating condensate in a harmonic trap) . 

The extreme rotation regime can be reached in two ways : either by letting tend to infinity with 



or by keeping il of the order of a constant (strictly larger than 1) and letting k tend to zero with 



where fio is a constant. 

In a recent series of papers [15l [161 El] > the large SI limit has been studied for traps related to the one we 
study (homogeneous traps, V{x) = \x\^ with s > 2, and flat trap V{x) — for < 1 and V{x) — +(X) for 
|a;| > 1). Their analysis applies to our problem in regimes "opposite" to the two described above, namely 
under the condition that either SI ^ G^^^k (for the regime opposite to (|1.44p ') or with k ^ G~^ (for 
the regime opposite to (|1.45p l. In this kind of regimes, the energy and density asymptotics are given to 
leading order by a limiting problem of Thomas-Fermi type. The limiting energy depends on the density 
only and is obtained from (|1.14p by neglecting the kinetic (first) term. 

If Slo is allowed to be 1 in (|1.45p . we have Dq, and so the potential term in the energy becomes 
negligible compared to the kinetic one. A solution should then converge to its projection on the first 
eigenspace of the Ginzburg-Landau operator — (V — iojx ) . This space is called the lowest Landau level 
and has been widely used for the study of rapidly rotating Bose-Einstein condensates in harmonic traps 
(see [5] and references therein). For the case SIq = 1 we refer to fSj where we have studied a regime of 
this kind in the lowest Landau level. 

The main drawback of the method we present here is that it confirms the appearance of the giant 
vortex state but does not give the critical speed at which the transition would be expected to happen. 
Finding such a critical speed is an important issue for the condensed matter community (see for example 
[251 133j ). For simplicity and comparison with other results we will consider k as fixed in the following 
discussion. 

At rotation speeds much smaller than those we have considered, the prefered state is known [T7] to be 
a vortex lattice encircling a central low-density hole. The transition between this state and the giant 
vortex is not expected to happen in the regime Sl ~ G^/^. Indeed, the results in 17] (see also [M]) allow 
to conjecture that in the case of a flat trap {V{x) = for |a;| < 1 and V{x) — -foo for |a;| > 1) the 
transition should happen when Q, ^ G(logG')~^. Combined with the scalings in [T^ this suggests that 
for our trap (jl.Sp the vortices will disappear from the condensate when ~ G'^/^(logG')~^, before the 



(1.44) 



n^^io>\, k^o, k<s:G-^ 



(1.45) 
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regime we have studied here is reached. 

In a forthcoming work [TH] we study, in the case of a flat trap, the transition regime H. ~ G(logG)~^ 
(corresponding to G^/^(logG)^^ for the trap p.Sp ) using tools from the Ginzburg-Landau theory 
[TTl HT] . This completes the analysis in [T7] where it is proved that the vorticity is uniform in the bulk for 
logG <C <C G(logG)~^ and bridges between the situations considered in [151 [HI [T^ and the present 
paper. We do believe that statements such as Theorem 11.21 and 11.31 are valid in the whole regime p.25p . 
but the proofs should rely on the tools used in [18J that allow for a better understanding of the energetic 
cost of the vortices. 



We make one last remark before turning to the proofs of our results, regarding the validity of the 
Gross-Pitaevskii description. It is expected to be valid when the number of particles is much larger than 
the number of occupied states. In that case a significant fraction of these particles must occupy the 
same energy state, which is the phenomenon of Bose-Einstein condensation. In Section 2 we analyze in 
some detail the eigenstates of the single-particle Hamiltonian (corresponding to the linear part in p.l6p ). 
In particular, we show fCorollarv l2.ip that the energy splitting between two neighboring single particle 
states is bounded below by a constant. On the other hand, the interaction energy per particle (remark 
that we have taken a unit mass constraint, so that actually represents the energy per particle of 
the condensate, and G /jja Im^^I^ the interaction energy per particle) is of order Gw^/^ as a consequence 
of Theorem 11.11 The number of occupied states can be computed by dividing the interaction energy 
by the energy splitting between two neighboring single particle states, so that one should expect the 
Gross-Pitaevskii description to be valid when 



where N is the number of particles. 

This is interesting because the Gross-Pitaevskii description for rotating bosons is known to break down 
in some regimes of fast rotation, even at zero temperature. It is the case for a condensate trapped by an 
harmonic potential of the kind (|1.3p when the rotation speed is too close (in some sense related to the 
number of particles considered) to the limit value min(ai, a2). In that situation, one must go back to the 
Hamiltonian for N bosons and the system exhibits strongly correlated states (fractional quantum Hall 
effect). We refer to [H^ and references therein for details on this phenomenon. 

The fractional quantum Hall effect happens in a regime where the minimization is restricted to the lowest 
Landau level. The corresponding situation in our setting would be a regime like p.45p but with llo = 1, 
which we do not allow in this paper. We are thus not in a regime where one should expect to get strongly 
correlated states. 



The article is organized as follows: in Section 2 we first gather the main notation that is to be used 
in the paper, then we study the one-dimensional energies (jl.lSp . providing the lower bounds on their 
ground states and first excited levels that we use in our analysis. In Section 3 we show how the results 
of Section 2 allow to bound the quartic term in (jl.lip from below and prove Theorem ll.il The proofs of 
Theorems 11.21 and 1 1 . 31 are then presented in Sections 4, with each main step occupying a subsection. 

2 The one dimensional energies 

We recall the definition of the one dimensional energies we are interested in: 




(2.1) 
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for every / G H'^ {R+ , rdr) n L'^(R+ ,rVn{r)dr). The potential F„ is given by: 



n 



Vnir) = y(r)+ --Lvr (2.2) 
r 



2nuj + (1 - Dn)Lo^r'' + y + y''^- 



2 



Where 



2 



y(r)=i?oy (r^-l)' (2.3) 
is the potential in the energy p.l4p . 

A straigthforward computation shows that Vn has a unique minimum point for r > 0, which we shall 
note Rn- It is uniquely defined by the equation 

Rl + ^-J^K = if-. (2.4) 

and it is interesting to note that ii n = lj then i?„ = 1, which is the minimum point of the original 
potential V . 

It is classical to show that there is an increasing sequence {^j,n) of eigenvalues for the operator 
associated to Fn- The corresponding sequence of normalized eigenfunctions {gj,n)j^i...j^^ is a Hilbert 
basis for the space H^iW^ ^rdr) n L^(]R+, rV„(r)(ir). Considering the one-dimensional functions gj,n as 
radial functions defined on , one has 

- ^9],n + Vngj^n = ^j,n9],n (2.5) 

and 

^ glJx)dx ^27: [ glJr)rdr = 1. (2.6) 
Jr+ 

We want to study the asymptotics of the first two eigenvalues Ai^„ and X.^n of the operator associated 
to Fn- We remark that (71 is unique up to a multiplicative constant of modulus 1, that we fix so as to 
have 

gi,n > 0. (2.7) 



Ai „ satisfies 



Al,r! — Fn{gi^n) 



and A2 . 



inf^F„(/),2^ / f{r)rdr^l\ (2.8) 



A2,n — Fn{g2,n) 

= inf |f„(/), 27r / f{r)rdr = 1, / f{r)gi^n{r)rdr ^O}. (2.9) 



Note that 52, n is not uniquely defined a priori. 

In fact it is not necessary to carry out the analysis for every n. Indeed, we have the following result 

Lemma 2.1. Let Ai^„ be defined by Ii2.8]) . There are constants a,b> and c > -s/G so that for 

\n — uj\ > a u?'^'^ 
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or equivalently 

\R„ - 1| > 

we have 

Ai,„>cw (2.10) 

Proof. For |i?„ — 1| > ciuj^^/^, one has 

Ai,„ >2tt [ Vn{r)\gi,nir)\^rdr > 2nVn{Rn) f |5i,nWprdr = K(i?rO > V{Rn) > C2L0 
for two constants ci , C2 . 

On the other hand one can see from formula (|2.4p that oj \n — uj\\s equivalent to 1 ^ |i?„ — 1|. For such 
values of n and i?„ we have — 1| ^ Cuj~^\rL — a;|, for a certain constant C, thus a; ^ |7i — a;| > 030;^/^ 
is equivalent to 1 ^ — 1| > C2(x'^^/^. We then note that the function mapping i?„ to n is strictly 
increasing and has a strictly increasing inverse, which allows to deduce that |rt — a;| > €30;^/^ is equivalent 
to \Rn — 1| > 02^"^^^ for any value of n and i?„. There only remains to tune the constants ci, C2, C3 to 
conclude the proof. □ 



AAi/2 = |neZ, \n-uj\<auj^'^^ (2.11) 



We will denote 



and study in detail the ground state and first excited level of the energy (|2.1[) under the mass constraint 
(12. 6p for n G A/1/2 • For n ^ the simple lower bound (|2.10p matched with the upper bound (|1.4ip 

will be enough for our purpose (note that (2Z?q + 4)^''^ < v^)- 

We now describe the main results concerning the modes n e A/1/2 • As we will see Qi^n 

and g2.n are 

asymptotically confined on a domain centered on i?„ which size is of the order of uo^^/"^ . More precisely 
we denote 

2 '-.^i' 



From (|2.2p and (j2.4p it is straigthforward to see that there is constant h so that for any n E A/1/2 

/i„ = /itj-i/2(i^_o(-l))^ (2.13) 
We define the rescaled functions ^i_„, i G {1,2} by 

(i,nix) = C~lgi^n [Rn + Kx) (2-14) 

where = 1 and 

2 _ 1 1 

Remark 1. In order for the asymptotics in equations |J.g7[ ) and h2.15\) to he justified, we need that ^ 
o-nd ^i^n decay faster than any polynomial when x goes to infinity. A typical example is to take for ^ a 
Gaussian, or more generally some function with an exponential decay. Very often in the rest of the paper 
we will have to deal with quantities of the form 

F{x)dx 
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with : R — !■ R an exponentially decreasing function (F[x) — x^e , fc G N for example), and 
dt^ = 0{lo'^), k > 0. We always replace such quantities by the corresponding 

/ + 00 
F{x)dx 
'OO 

which is justified because the decay rate of F guarantees that the difference between the two will be expo- 
nentially small as to goes to infinity, whereas all the other remainder terms that we will encounter will 
behave like powers of lo. 

Note that the decay o/^i.„ is proved in Proposition li2.S^) . 

Notation. In the rest of the paper, we wiU always mean by 0(w'^) a quantity bounded by w'^ uni- 
formly with respect to Dq and n. C will be a generic positive constant depending neither on lo or Dq 
nor on n. When writing that a function is a 0^2(0;'^) for example we mean that the function's norm 
is bounded by w*^ uniformly with respect to Dq and n. For the rescaled functions of the type (|2.14p 
we will often use the notation Oho{^''), the HO symbol meaning that the 0{u)^) is taken in the norm 
associated with the harmonic oscillator p.24p . namely the norm in H^{R) n L^(R,x^dx). 
We recall that the notation a oc b has the usual meaning that a/b converges to some constant. 

From now on and in the rest of this section, we always implicitly consider only those n that belong 
to A/1/21 which have the following important properties deduced from (|2.2p and (|2.4p : 

• In — Lu\ < aLU^/"^ 

• Rn and are bounded uniformly with respect to uj and Dq, more precisely — 1| < buj^^^'^. 
Moreover we have 

\Rp - Rq\ CX uj^^\p-q\ for any p,qe 7Vi/2. (2.16) 

• hn OC w"^/^. Moreover we have 

\hp - hq\ CX - q\ for any p,q £ Mij2- (2-17) 

• T/„^'^'(i?„) CX up' for any k >2. 



We are able, using the localization property of f;i_„ and 32, n, to expand the corresponding eigenvalues 
in powers of /i„. We have the foUowings results: 



Theorem 2.1 (Asymptotics for the ground states of the one-dimensional linear problems). 
Let Xi^n be defined by i2.8\) . 

Suppose that n G A/1/2 ^5 defined in equation \2.11]) and let Rn be defined in equation {2.4^ . 
We have as lu 00 and Dq — > D 



\l,n = Vn{Rn) + \] "'^ +Kn (2.18) 

where Kn = 0(1). Let ^1 be the normalized ground state of the harmonic oscillator {L24^ and ^i_„ be 
defined by equation \2.14^ . We have 

ii,n = a + KPn^i + hlQnii + OHo{io-'"^) (2.19) 
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where /i„ is defined in equation \2.12\) . Pn and Qn are two polynomials whose coefficients are bounded 
independently of uj and Dq and depend continuously on n in the following sense: 

Pp^i = P,^i+OHoito-'\p^q\) (2.20) 
Qp^i = QqCi+OHoi^~^\p-q\) 

for any p,qG J\fi/2- 
Note that 

K.(i?„) + /^^ = 0(-). 
From this theorem we deduce as a corollary 
Corollary 2.1 (Variations of Ai.„ with respect to n). 

Under the assumptions of Theorem \2.1\ and for uj large enough, Ai.„ is minimized for at most two integers 
n* — Lu + 0(1) and n* + 1. For any n ^ n* , n* + 1 one has 

Ai,„ = Ai,„. + C (n - n*f (l + 0(c^"i/2)) . (2.21) 

Equation (j2.2ip points out the dependence of Ai_„ with respect to n, which is the key point to prove 
Theorem ll.il We now state the result giving the gap between A2,n and Ai,„: 

Proposition 2.1 (Lower bound on the first excited levels of the one-dimensional energies). 
Let A2,n be defined by \2.9\) . 

For every n e A/1/2, ™c have as u +00 and Dq — > D 



X2,n > K(i?„) + 3^^^- > ^ 2^1^^. 



(2.22) 



Note that it is feasible to expand ^i^„ further and to make of Proposition 12.11 a statement as precise 
as Theorem 12. 11 two improvements that we do not need in the sequel. 

For convenience we sum up some properties of gin that will be useful in the rest paper in the following 
Corollary. 

Corollary 2.2 {L^ bounds for (?i,n)- 

For any n G J\fi /2 and any 1 < p < +00 there is a constant Cp depending only on p such that 

\\9l,n\\LP(R+,rdr) <GpUj'^'^^^ (2.23) 

with the convention that ^ — if p — +00. 

Proof. These estimates are simple consequences of (j2.14L (j2.15p and (j2.19p . The integrals are computed 
using the change of variables r = Rn ~\- hnX. Remark that ci^n < Cuj^/^ with C independent of w, Dq 
and n. □ 

We proceed to the proof of Theorem 12.11 assuming that its assumptions hold true for the rest of the 
section. We begin with an upper bound on Ai^„: 

Lemma 2.2 (Upper bound for Ai_„). 

Under the assumptions of Theorem \2.1\ we have 



Kn < K(i?„) + /^^+ X ^'^'^^^^^ + ^'-'^^ 
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Proof. We proceed as for the proof of (ll.4ip . taking as a trial function 

9T\r)^c^i{h-'ir-R^)) (2.25) 

n \\L^{m+,rdr) 

1/4 

(2.26) 



where c is chosen so that 27r|| (7^'^^*||?2|-»+ = 1 with 



An important fact for the calculation is that for the range of n we are considering we always have 
hn = 0{uj~-'-/^). We also use that is an even function. 

□ 

We now show the GxpoiiGiitia.! decay of yi,n and its derivative outside of a shrinking region: 

Proposition 2.2 (Pointwise estimates for .gi,„). 
Let gi^n be defined by \2.8\) . 

There are positive constants a\ and ai so that for every n G A/1/2, every r G and for every e > 

|ffi,„(r)| < C,uj^'^+'e~<'^y (2.27) 

and 2 

|g'i^„(r)| < CeC.3/4+ee-=(^) . (2.28) 

Moreover, for any a < 1 there is a constant Ca so that 

\9iAt)\ < Ce-^'^r^ (2.29) 

for any r < Ca ■ 

Proof. The proof is done in four steps. 

Step 1. We first need a global L°° bound for .gi,,i. We claim that 

Vr/>0, Ve > 0, 3Ce > so that ||gi,„||L=o(],,^+^[) < Cewl/4+^ (2.30) 
We first note that ipTM]) implies 



Al,n = i^„(gi,„) < Vn{Rn) + (Rn) (2.31) 

because both T/^'(i?„) cx and V^^\Rn) oc when n G A/1/2 • Thus, using (|2.3ip and (|2.6p : 



2^ / (l5'i,„WI' + (K(r) - VniRn)) \gi.n{r)n rdr < VWiRn 
Jr+ 

but Vn{r) -Vn{Rn)>0 SO 

2TTrj f g'ljrfdr < 27t [ \glnir)\''rdr 
< 27: [ \g[Jr)\''rdr 
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and recalling that F^'(i?„) = 0(t^^) for n G A/1/2, using again 

\\gi,n\\HH]n^+oo[)<Cuj^'\ (2.32) 
Now, an interpolation argument, using (|2.6p and (|2.32p yields 

Ve>0 ||gi,„|L 1^3.^, ^ <CeC.i/4+e (2.33) 

H 2 (]i;,+oo[) 

which implies (j2.30p by Sobolev imbedding. 
Step 2. We have 



and Vn{Rn) + ^V^'(i?n) = 0(cj), thus there is a 7 > so that 



Vn{r) - Vn{Rn) - VW{Rn) > for |r - 1| > 7c^-i/^ (2.34) 
Using that i?„ = 1 + 0(w^^/^) and /i„ cx oj^^/^ we deduce that there is a (5 > so that 



Vnir) - V„{Rn) - VVniRn) > for \r - i?„| > SK. (2.35) 

We denote 

Sn = {reR+, \r-Rn\<SK} (2.36) 

and prove an exponential decay property for gi^n on 5^^. We recall the Euler-Lag range equation for (71, n 
(considered as a radial function defined on R^) 

- Agi^n + Vngiji = Ai,„gi,„. (2.37) 

Using (|2.3ip we have for lo large enough 

- Agi,„ + {Vn{r) - VniRn) - ^V^iRn)) ffl.n < 0. (2.38) 

We now apply the comparison principle on 5^ : the right-hand side of (j2.27p is a supersolution for (|2.38p 
on 5*^ with a proper choice of cji, the boundary condition being fulfilled thanks to (|2.30p . Thus the result 
holds on 5^. On the other hand, on Sn the inequality is true for large enough. Indeed, we have the 

global upper bound (|2.30p and on 5„ the function e is bounded below by some constant. 

Step 3. We prove (|2.29p . It is straigthforward to see that the function on the right-hand side of (|2.29p 
is a supersolution for equation (|2.38p on Bc^ for some constant Ca depending only on a. We then use 
the maximum principle on Bmin(CQ J-) f*^ conclude that (|2.29p holds. The boundary condition is fulfilled 
thanks to ^m\ . 

Step 4- We now prove (|2.28p by an interpolation argument. Using (|2.37p . (|2.24p . the definition of Vn, 
f^n7\i and l^r^i we have 

|A5i,n(^)| < Ccj^/^+"e-"-^(^^)' (2.39) 
for some S3. Using Gagliardo-Nirenberg's inequality (see for example pUl HO] ) 

||Vgi,„||L~ < C||Agi,„||^(^||5i,„||^{^ 
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on Br for any < i? < i?„ and on i?^ for any i? > i?„, we get from (I2.27P and (|2.39p 

for some S2, which imphes (|2.28p . 

We recall the definition of the blow-up function ^i.„: 

9i.n{r) = ci,„^i,„ {h~^ {r - i?„)) 
where /i„ is defined in equation (|2.12p . 

Note that this function is only defined on ] — +00 [, but as —j^ —00 and gi,„ decreases exponentially 
fast in this limit, we may (remembering Remark 1) abuse notation and consider it as a function defined 
on M. 

The next lemma states a few estimates that will prove useful to expand Ai.„. They are consequences of 
Proposition [ 



□ 



(2.40) 



Lemma 2.3 (Estimates for 

Let n be in A/1/2 fl*^*^ C'fe.s generic constant depending only on k and e. We have, for lo large enough 
and every e > 

^- |ci,n^n J^x''£,i^,i{x)^ {Rn + hnX)dx\ < Ck^e^" ■ 

3. \hn!^x''^i,n{.xfdx\ <Cfc,eW^-l/2. 

4. \h^j^x^^'l^{x)dx\<C^^,u^-^/\ 

Proof. We start with the proof of 1. Using (|2.27p we have for every e > 
Ci„/l„ / x^£,i^„{xf {Rn + hnx)dx : 



r- R, 

+ V hn 



glA^ydr 





r-Rn 


/ 


hn 



'/^+'e-^''i'^)\dr. 



A change of variables r = i?„ + hnX then yields 

fe 



^ Rn \ 2 



/R+ V hn 

We turn to the proof of 2. We have 



9i.n{r)rdr 



< C,Lu^ I x'^e-^'^'' dx < Cfcecj^ 



'^l.nhn — Ci nhn / 



R+ 



9i.n{r)dr 



(2.41) 



(2.42) 
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but 



Rn 



9i,n{r)rdr + — 

R+ -ftri 



{Rn -r) glJr)dr 



1 1 



+ 15- / (Rn ~ r) gl„{r)dr 
{Rn-r) gl„ir)dr. 



_^ 1 

Rn 



Using (|2.27p again one gets 



Now 



{Rn - r) glJr)dr 



hn gi^nir)dr < {Rn - r) g-^^^{r)dr < K l gi,„(r)dr 



(2.43) 
(2.44) 



and gathering (|2.43p and (|2.44p we get 



9l,n{f)dr- 



2TTRn 



for some constant C. Using (|2.42p we have proved 2. 
The proof of 3 uses the same kind of computations: 



/ x''£,ln{x)dx 



1 f fr -Rn 



9i,nir)dr 



< Ck,6^'^hn. 

We have used (I2.27P for the first inequaUty and the point 2 for the last one. We prove 4 exactly like 3, 
using (p:^ instead of (PTTT]) . 

□ 

We are now able to present the 



Proof of Theorem \2.1l We split the proof in four steps. 

Step 1. We use the blow-up (|2.40p and expand Ki(r) around i?„ in the expression of Ai^„. Using the fact 
that 

VjeN, VneAAi/2 V'^^\Rn) ^ 0{u^) (2.45) 

and hn = 0(cj~^/^) we get 



Ai,„ = F„(5ri,„) = 27r-^ / £,[ „^{xf {Rn + hnx)dx 



hn 



fc=0 



h'^+^x''^l,n{x)^ {Rn + Kx) dx 

+ 0(c^2-Ar/2) /■ UA^^^ „(a;)2 (^^ + /i„x) I (2.46) 
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where 



'l,n 



27r/l„ / Cl,nix)^ (Rn + hnX) dx 2TThnRn + 27r/l2 X^i^nixYdx 



(2.47) 



Then we use the estimates of Lemma [^751 and (I2.45P again to bound the coefficients of Ai^„ in its expansion 
in powers of /i„ (|2.46p . Recahing the definition of /i„ (|2.12p we obtain 

Al,„ = Vn{Rn) + \f^^l^ i^lnixf + X%A^f) dx + 0{u^^^+'). (2.48) 

for every e > 0. Combining this estimate with (|2.24p and dividing by y^-"^^^^ = 0(w) we ffiid 



^ iCii^r + ^'u^r) + o{u-^'^+^) > ^ ieu^r + x^^.^^xr) > £ {e,{x)' + x^^,{xr) , (2.49) 

but ^1 is the unique normahzed minimizer of the energy associated to (|1.24p . so we deduce 

^i,n=^i+OHo{uJ-'^^+n (2.50) 

for every s > 0. The next steps consist in improving (|2.50p . using first a sharper energy expansion then 
an equation satisfied by 

Step 2. We progressively improve the remainder term in (I2.48P by a bootstrap argument. Going back 
to (|2.46p and using again the results of Lemma [^31 we may be more exphcit: 



Vr'/iRn) 



Xl,n > VniRn) + ]J ^^^^ (^L„(^)' + ^^^lA^f) dx 

+ Clu^^^+'[ [ x^i^„{x)^dx+ [ x'^ii.r.ixfdx 
for any e > 0. If, for some /3 we have 

a.n=Cl+OHo(^^), 



x^[,^{xfdx]-Cuj' (2.51) 



(2.52) 



then it is easy to show that for every £ > 



x^^iAx)^dx - / x^^i{x)^dx 



(2.53) 



using (|2.52p . (|2.27p and the fact that is a Gaussian. Now, being an even function we deduce 



a;^Ci,n(a;)^dx 



(2.54) 



and arguing hkewise to bound \J^x£^i^n{x)'^dx\ and „(a;)^da;| we improve (j2.5ip into 



> + \f^^l i^lni^)' + ^'a.n(^)') dx + C^'/'+^+f (2.55) 



and the same argument as that used with (|2.49p to get (|2.50p yields 



(2.56) 
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1 P e 
'4+2 + 2 



for every e > 0. The fixed point of tlie function 

being /3 = — ^ + e, we deduce by induction that for every e > 
using (j2.50p as a starting point for the induction. 



(2.57) 



Step 3. We are going to improve the expansion (|2.57p using the equation satisfied by We claim 



that 



^l,n - a + KPnil + Oho{0J-'+') 



(2.58) 



where P„ is an odd polynomial of degree 3. Let us write the equation for a,„. We write (j2.37p in radial 
coordinates: ^ 

- - + K Wffi,nM = Ai,„gi,„(r). (2.59) 

Then, making the change of variables r ~ i?„ + hnX, expanding Vn and ^ and multiplying by h'^Ci]-^ we 
get 



Ai,„/i^,a,n + 0L2(Lc>-i+^). (2.60) 



We have used the estimates of Lemma [2.31 again to evaluate the remainder. Taking /3 = — 1/2 + e in 
(|2.55p we have 



so, using the definition of /i„ ((2T2)) . ((2^ reduces to 



S "n^ = + )■ 



Note that fe^ = 1 by definition. Now, if we write as 

6,n = + ^n<y5n 

and insert this in (|2.6ip . using the Euler-Lagrange equation for and the estimate 

(Pn = OHo{^') Ve>0, 
deduced from (j2.57p and (j2.13p . we get an equation satisfied by (p„, namely 

On the other hand, the mass constraint for ^i.„ allows us to deduce from (I2.62p 



(2.61) 



(2.62) 



(2.63) 



(2.64) 
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Using an expansion on the basis of normalized eigenfunctions of the harmonic oscillator (|1.24p . one can 
see that the problem 

r -u" +,^u-u^^j[^ ^^hix^i. 



has a unique solution in iJ^(R) n L^(M,a;^dx), which is of the form u = Pn^i where P„ is an odd 
polynomial of degree 3. Substracting (|2.63|) from the first line of (|2.65|) and using (|2.64[) it is easy to show 
that 

= Pn^l + Oho{^'^'^^') (2.66) 

thus proving (|2.58p . The first line of (|2.20p follows by inspection of the right-hand side of (|2.63p . 



Step 4- We complete the proof of (|2.19p . The argument is similar to that of Step 3. We can compute 
from ((^351) that 

Xl.n^Vn{Rn) + ^f^^+Kr,, (2.67) 

where Kn = 0(1). We write 

il.n = a + hnPnCl + /l' V'n (2.68) 

and arguing as in Step 3 we get an equation for ipn'. 

+ - ^^^^Kx^P^^, + Oi.(a,-i/2+e) (2.69) 

and the condition 

j eiV-n = -\j P?A! + 0{u-^'^). (2.70) 
Note that, multiplying (|2.69p by integrating and using (|2.70p we get 

if„ = < + 0{uj-^'^+') (2.71) 

where 

^ ^ Ri^^ ' 4! '""■^ i?„ ' 4! 

Solving ((2J9)) we find 

= Qn^l + Oho{^~^'^^') (2.73) 

for every e > 0, where Qn is a polynomial and QnS,i is the unique solution of the problem 

This concludes the proof of (|2.19p . We deduce (|2.18p by straightforward calculations. □ 
We present the 



K fii( -^e, + - + ^^^^hix'P^^,] . (2.72) 



(2.74) 
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Proof of Corollary \2.1\ At this stage, we know that 

Ai,„ = C{Rn) + if„ 

where the function C is given by 



C(i?„)-K(i?n) + y^^^^^^ (2.75) 

and Kn = 0{1). Note that Kn — Km = 0(w^^/^) for any m,n € A/1/2, as a consequence of (I2.20p . which 
allows to focus on C. It is a function of i?„ only, remembering the one-to-one relation (|2.4p between i?„ 
and n. Recalling (|2.2p we have explicitly: 



C{Rn) = Dn'^Ri + (1 - t^'i?' + DnLo^ - 2u^^,RIJrI ' ^ 



+ ^12Dnuj^Rl + 12Lj^il- Dn). (2.76) 

Studying the variations of C(r) with respect to r G we can see that C is minimized at a unique 
R S R"*" depending continuously on Dq and We have no explicit expression for R but one can see that 
R= 1 + 0{uj^^). Expanding C around R we have for any n g A/1/2 

C(i?„) = C{R) + (i?„ - i?)' f 1 + Oiiu-'/')) (2.77) 



2 

because C^'^^R) = 0(a;^) for k > 2 and i?„ — i? = 0(cj^^/^) when ?i e A/'i/2. On the other hand, using 
(EH we get that for any n G 7Vi/2, w|i?„ - i?| cx - n| where N = dI^^uj (^i?^ _^ i^i?^^ " 

C(i?„) - C{R) cx (n - Nf (1 + 0(c.-i/2)) . (2.78) 

We want to define n* as the minimizer in A/1/2 of ("^ ~ A^)^. Only two cases can occur : either A'^ is 
exactly half an integer, then [A^] and [A^] + 1 both minimize (n — A^)^ in A/1/2, or there is exactly one 
minimizer. 

In both cases we pick one minimizer n* and (I2.2ip follows by inspection of (j2.78p . An important point is 
that \n* -N\ < 1/2. 

□ 

We now turn to the proof of Proposition l2.1l We begin with an upper bound on A2.„: 

Lemma 2.4 (Upper bound for A2.„). 

Let A2.n and 52, n be defined by i2.9\) . 

We have for every n G A/1/2, as co ^ 00 and Dq D 



(2.79) 



A2.. = F„(g2.„) < K(i?n) + 3^l^^+0{uj'^'). 
Proof. We begin with a test function constructed as that in the proof of Lemma 12.21 

GT'{r)^c^,{h-\r^Rr,)) (2.80) 



where c is chosen so that 27r||G^'^**||^2 = 1- Then we define 



= GT - / giAr)GT\r)rdr gi,„ (2.81) 
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and, as /j^^ G^^^* gi^nrdr = we have 

Fn{92,.) < F^r). (2.82) 

Using the Euler-Lagrange equation (|2.37p for gi^n one proves that 

F„(Gr*) = F,,{GT') ~( [ g,^r.{r)GT'{r)rdr] Ai,„ < F,,{GT')- (2.83) 



R+ 



Computing F„(G^^*') as (|1.4ip . we obtain 



= Vn(R,,) + 3^l^^+0{oj'/'). (2.84) 

because ^2 is the second normalized eigenfunction of the harmonic osciUator p.24p . associated to the 
eigenvalue 3. On the other hand 



^ 2 

^test I 



\\GT\\U^^,r,r) = 1 - [J^^ 9iAr)GT\r)rdT ^ 

>l-c(^j ^i.n{x)£,2{x){Rn + Kx)dx^ >l-Cuo-^ (2.85) 
using (j2.19p and /]g^iC2 — 0. Gathering equations (|2.82p to (|2.85p we get the result, remembering that 

□ 

Proof of Provosition \2.1\ Starting from the upper bound (j2.79p . we can prove the equivalents of Propo- 
sition [511] and Lemma for g2,n- We omit the detailed calculations since they are easy modifications of 
those we used for gi^n- Just note that 32, n is not positive so we have to rely on the inequality 

- A|52,„|' + 2Vn\g2,n? < 2A2,„|g2.n|' (2.86) 



to prove the equivalents of (|2.27p and (|2.29p . (|2.86p is a consequence of (j2.5 
Then we can argue as in Step 1 of the proof of Theorem 12.11 and obtain 



A2,„ = VniRn) + \f^^j^ (?2,„(^)' + x'^2Axf) dx + 0{lO^'^+'). (2.87) 

But 

/ gi,n{r)g2,n{r)rdr = 

and this implies, using (j2.19p . 

' ii{x)i2Ax)dx = 0{uj-^'^) (2.88) 



We obtain from (PTST)) and ^I^. The second inequality holds because Vn{Rn) =0{uj'^). □ 
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3 Lower bound on the interaction energy of m^, proof of Theo- 
rem 11.11 

In this section we show how the estimates of Section 2 on the one-dimensional problems allow one to 
bound the quartic term in the energy p. lip from below. We then deduce the energy and density asymp- 
totics of Theorem ll.il 

We begin with a property of exponential decay for u^: 

Proposition 3.1 (Exponential decay for u^). 
There is a constant a so that 

\u^{x)\ < Ct^i/^e-'^'^ll^l-il' for any x G (3.1) 
Proof. We recall the Euler-Lagrange equation for Uj^: 

2 

- (V - iujx-^fu^ + - 1)^ u„ + 2G\u^\'^u^ = fi^u^ (3.2) 

where /i^j is the Lagrange multiplier associated with the mass constraint. From this we deduce an equation 
satisfied by U := \u^\'^, 

-AU + AGU^ + 2 |Vu^ - iiux^u^l^ + 2U (^i'o y (NP - l)^ - fij^ - 0. (3.3) 

As in the proof of Proposition 12. 21 we need first a global L°° bound. We claim that 

\\U\\l^ < Coj, (3.4) 
which is proved using a method due to Alberto Farina [20j . We introduce 



^-^-2G- 



We use Kato's inequality and 



A(C/+) > %>oAC^ 



= l^>o {^GU^ + 2 \Vu^ - iu:x^uj^ + 2[/ (^^o y (l^^l' - l)' - A 

> %>o (4GC/2 - 2^^[/) 
= ly>o4GC/2 > 4G(C/+)^ 

We have -A(C/+) + 4G(C/+)2 < which implies [7+ = (see [13]). We multiply the Euler-Lagrange 
equation (|1.16p by u*^ and integrate over to get 

M„ = /^ + G / \u^t (3.5) 

and thus 

Ma. < < Gw (3.6) 

because of the upper bound (jl.4ip . Using (j3.6p and the fact that G is larger than some constant, the 
claim is proved. 

The rest of the proof is an application of the maximum principle. We define 



\x e M^ - \f < Suj-^^ . (3.7) 



24 



We have 



,|2 



^U + 2U[Dn—{\xY-ly-^i^]<Q (3.8) 



and (^D^^^ (|a;p — l)^ — ^ij^ > on A^^^^ if 5 is large enough. We use two comparison functions: 

which is a supersolution to (|3.8p on a H < 1} if ci is chosen small enough and 

which is a supersolution to p.8p on A^^ 5n{|a;| > 1} if cr2 is chosen small enough. The boundary conditions 
for the application of the comparison principle are fulfilled thanks to (|3.4p by taking Ci and C2 large 
enough. There remains to note that 

Fi{x) < Cic^e-"i"(l"l~i)' 
to show that (|3.ip holds on ^ with cr = ^ min((Ti, (72) and C = ma.x{^/Ci, \fC^. 

On the other hand (j3.4p implies that (j3.ip is true (for a large enough constant C) on Aujfi because 
^-au){\x\-i) bounded below there, which concludes the proof. □ 

Let us recall that we have the Fourier expansion 

We denote (., .) the scalar product in L^(M+, rdr): 

{u,v) = / u{r)v*{r)rdr 
Jr+ 

and introduce 

u^(r,0):= (/«-5i.n>5i,n(r-)e™^ (3.9) 

n6AAi/2 

As a consequence of our analysis in Section 2, the modes corresponding to n ^ A/1/2 carry very little 
mass, whereas for n € A/1/2 we have /„ cx gi^„ in the sense. This can be summed up in the 

Proposition 3.2 (First estimate on \\u^ — Uu)\\l^{r2)). 

Let Uuj and u^j he defined respectively by \1.15\) and h3.9(l . Under the assumptions of Theorem \1.1\ we 
have 

11"^ - u^Wmv?) < CGi/^u;-!/^ = o(l). (3.10) 

Proof. We have to prove 

V / Ifnl^rdr < CGlu-'^' (3.11) 



and 
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For n G A/1/2 can write /„ as 

1/2 



f,i = {fn,gi,,i)gi,n+ (^J ^\fn~ {fn,gi,,i)gi,n\'^rdrj (3.13) 
/ gi,n{'>')gn{r)rc 

and 

2^ / Is^MP'^rf'^ - 1, 



where 

rdr = 



thus we have 



-Fn(/n) > Ai.„ |(/„,5i.„)|^ + A2.„ / l/n - (/n,5l.n)ffl,n|^rdr (3-14) 

for n e 7\/i/2. Using (|1.43p and p.l4p . ignoring the interaction energy, we bound F^{Uijj) from below in 
the following way 

Fu,{uu,)> ^ Ai,„ |(/„,5i,„)|^ + ^ A2,„ / \fn- {fn,gi,-n)9iA^ 

Ai,„27r / \fn\^rdr. 



But the energy of gi,„»e™ ^ is an upper bound for Fi^[uui) and Ai_n > Ai_„* for any n. Combined with 
the results of Proposition 12.11 and (|2.10p we get 



Ai,„. +27rG/ gtn'ir)rdr > V (|(/„, gi,„)|' A 
JR+ „f77 ^ 



+ |^Ai,„- + Y ^^-^ j y + ~ (/n,ffl,n)5l,n| ''drj 

+ V cw27r /" |/„|Vdr (3.15) 
with c > \/6- On the other hand, we know from the normalization of u^^ that 

neZ '' neJ\fi/2 

+ X! \fn- {fn,gi,n)gi,n\^rdr 

+ V 27r / |/„|Vdr. (3.16) 
Combining p.lSp and (|3.16p . using that Ai_„* — ^/6uj < for w large enough we get 



27TG f gl^,{r)rdr> ^ / | /„_(/„, g,,^^,,^^ ,/ ^^^^ + ^ Cc.27r/ |/„|Vdr 
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and this implies the desired result because \/ ^" o^"'' > Cuj for rt e A/1/2 and 



which is a consequence of (|2.19p . 



□ 



We now prove a simple lemma which states that the total interaction energy of is an upper bound 
to the pairwise interactions of the modes /„. We state the lemma for a general $ € L^{M?): 

Lemma 3.1 (Generic lower bound on \^\'^)- 

Let $ be any L*(R^) function with Fourier decomposition 



We have 



Proof. We expand |<I>|^: 



$(r,0) = ^0„(r)e"^ 



(3.17) 



m 



E 



+ Jo 



(3.18) 



and the only integrals with respect to 6 that are not zero are those for which n+p~q — r — Quo 



I \<^\'^ = E / 4>n+q(f>p4'q4'n+prdr 



/ 



> 2tt 



rdr 



S^rV / \<j>p\^\cj)g\'rdr 



□ 



We are now able to prove the lower bound on J^a combining the results of Proposition 13.21 and 

Lemma 13.11 A very important fact in the computation is that for n e A/1/2, fn oc gi.m and that the 
appropriate blow-up of .gi.„ converges to a gaussian (equation (|2.19p l. 

Proposition 3.3 (Lower bound on the interaction energy). 
Let Uuj be a solution to hl.l5\) with Fourier expansion 



E/"W' 



,i7i0 



Set 
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where gi^n is defined by equation We have 



+ 0{Lu'^^\\u^~ujL-) + Oil). (3.19) 

Proof. We begin by using Lemma |3. II 

/ \u^\'>2nJ2 I \fv?\U?rdr. 

A consequence of Proposition 13.21 is that it is enough to consider only the pairwise interactions of the 
modes labeled by n S A/1/2 (see equation p.38p below), so that we use 



and evaluate the right-hand side. 
We note that, using 



< 



< 



C(«(AAi/2))'e-^'^ < e^^". (3.20) 

We have used that tl(A/i/2) < Cuj^^"^. A similar estimate holds true for J2pqeJ\ri/2 fr<i/2 l/pPl/?!^'"^^ 
that it is sufficient to consider integration domains where 1/2 < r < 3/2, which we implicitly do in the 
rest of the proof. 

It is convenient to introduce the following quantities: 

an = Rn - Rn' (3.21) 

and the rescaled functions (p defined for p e A/1/2 by 

fp{r) = dj,Cj, {h;\r - Rp)) (3.22) 
where hp is defined by equation l|2.12p . 

4 = , ^ '{^''^'^ (3.23) 

hpj\Cpix)\ {Rp + hpx)dx 

and IICplU^ = 1. In this proof we denote /?„ the n-th term of a generic sequence satisfying 

J2pI<C (3.24) 
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where C is indenpendcnt of lo and D^. The actual value of the quantity /3„ may thus change from line 
to line. 

Using (|3.10p . we have 

V \dpCp - ci^p {fp, gi,p) ^i^pf {Rp + hpx)dx = V / \fp - ifp^ 9i,p) 9i..pf rdr 



< Cuj^'^\\u^-u^\\l. 

thus 



Vp e 7V1/2 dpCp = Ci,p {fp,gi^p) 6,p + OL2((flp+/ipa;)rf2;)(/3p^^^^''||Maj " 'iia;||L2) (3.25) 

where '^pPp is bounded uniformly with respect to lo. 

On the other hand, using the bounds ^n{fn) < -fw(ua)) < Co; we have 

\\dpCp\\m({Rp+hpx)dx) < /3pW^/''. (3.26) 

Also 

\\ci,p(l,p\\HmRp + hpx)dx) < w^/'' (3.27) 

follows from (|TTg|) . 

A change of variables r = i?„. + hn'X yields (q;„ is defined by (|3.2ip ) 



(3.28) 



Then we have, using p.25p . (|3.26|) . p.27p and Sobolev imbeddings 



^\{fp,9l,p)\^\{fq.9l.q)? K'cl^pCl^q j ^Ip (j^^+jf-^^lq {jT^^It) +^"'^)'^^ 

Uuj — Uuj\\l^)- (3.29) 

Remark that it is clear from the definition of ci^„ (|2.15p and (|2.19p that 

^^--2^ + ^(-"'') ^'-''^ 
because is an even function. Using (|2.19p and p.30p we obtain 

+ 0(a;i/2||y^_ (3 3^) 
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But for any p £ 7Vi/2 we have from ((TTTI that = 1 + 0{iLj~^\n* ~ p\) = 1 + 0{uj~'^/'^) because 
n* =LU + 0{l),so 



and using that is a Gaussian we have 



= (l - Ct^ {ap - aqf^ J £.t{x) (i?„. + hn-x) dx + ©(w^^/^) (3.33) 

where we have used that /i„ — 0(cj^^/^) and |a„| = |i?„ — i?„. | cx cj"^]?! — 7i*| = 0(lj~^/^) (see (|2.16p '). 
On the other hand we have 

Vp, q e -i- = -L + 0(^1/2) (3.34) 

using that ^ = 1 + 0(w-i|n* = 1 + 0{uj^^^^) and 

using that \Rp — i?„.| oc ijJ^^Ip — n*\ ~ 0(ti>^^/^). Gathering equations (|3.28p to p.35p we have 



1 



/ \fpW,\^rdr> \ifp^9i,p)\'\{f„gi.,] 

X (^-Cuj {ap - aqffj J ^^{x) {Rn- + K-x) dx 

+ 0{u:^'^\\u^-u^\\l2) + 0{1). (3.36) 

But ((2J9)) gives 

27r [ g,,r.'{r)^rdr= ^ f ^,{xf {R,,, + h,,,x.) dx + 0{uj-^^^) (3.37) 
We also have, as a consequence of p.lOp and (jl.lOp . 



2 

2 



(1 - - uu^h-f > 1 - C\\u^ - u^\\l2 (3.38) 
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and 

< 4 E (3.39) 
because EqeT^^/^ |(/„5i.9)l^ < 1- Combining djmi . ll^^ and 

we obtain 

•^K' pG.^i/2 

We conclude the proof by recalling that |q!„| — \Rn — i?n» | oc uj^^\n — n* \ using (|2.4p . so that 

□ 

Note that the n-th term of the second line of (|3.19p can be interpreted as an upper bound to the gain 
in interaction energy per unit mass carried by the mode /„. On the other hand Corollary 12.11 provides 
a lower bound to the loss in kinetic and potential energy per unit mass carried by the mode /„. The 
comparison of this two quantities is the crucial argument of the 

Proof of Theorem We proceed in three steps 

Step 1. We first claim that the following improvement of (|3.10p holds 



Ml^ <CGuj-^/^. (3.40) 
Combining the results of Theorem 12.11 and Proposition 13.31 we get 

FUu..) > 2nG [ \gi.r,.\''{r)rdr 

+ E l(/"-5i,n>l'(Ai,„-CGw-^/'|n-n*p) 

+ A2,„ / \fn~ {fn,gi,7i)gi.n\'^ rdr 

V Ai,„27r / \fn\^rdr + 0{Guj^/^\\u^-u^\\L2) + 0{G). (3.41) 



We now use (7i.n*e™ as a test function and proceed as in the proof of Proposition 13. 21 With the result 
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of C0r0llarv r2.ll we get 

C{Guj^'^\\u^ - u^\\l2 +G + Gu-^'^) 

> E \{fn,9i.,n)\^{c\n-n*\^-Gu;-'/^\n-n*\') 



E 



E C*^ / \fn - {fn,gi,n)gi,nf rdr 

y2 Clu [ |/„prdr. (3.42) 



Now, for n e A/1/2, one has 

C|n - n*|2 - Gu;'^^^\n -nf>0 
because G <^ uj^^"^, so equation p.42p yields 

G(Gco'/^\\u^-u^\\L2+G + Gio-'/A > Cuj y2 f \fn - {fn,gi,n)gi,n\^ rdr 



1/2 



L2 



which implies (|3.40p . 

Remark also that we can replace A2,„ |/„ - {fn, gi,n)gi,n\^ rdr by X]j>2'^i."l ifn^gj.n) P in p.4ip . 
Indeed, writing 

/« = E 93,n (3.43) 

one has 

Fnifu)^Y.\(f^^9j.n)\'Xj,n. (3.44) 

This yields 

G{G' + G + Glo-'/^)> E E^:'^"I(/"'5,.«>I' (3.45) 

which will be useful in Section 4. 

S'tep ^. We prove (|1.26p . The upper bound is obtained by taking gi^n-e™ ^ as a test function. For 
the lower bouind we first prove the following essential estimate: 

E \\fn\\l^{^+..rdr)W-n\^ <CG\ (3.46) 

We recall that we assume G^ > g^ for some constant g so it suffices to prove 

E \\fn\\h^u+,rdr)\n-n*\^<CG^. (3.47) 
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Coming back to (|3.42p and using p.40p we get 

c{g-' + G + Guj-'/^)> ^ |(/„,5i,n)|'(|n-n*|2-Gc.-i/2|n-n*|2). (3.48) 

But GiLj~^^^ ^ 1, so if 71 7^ n*, 71* + 1 we have, for some constant C < C, 

C\n - 7i*|2 - GLo-^/'^\n -~n\^> G'\n - n*f. 
This allows to deduce from (|3.48p that 

E \{fn,gi,n)\'\n^n*\' <CG\ (3.49) 

On the other hand we already have as a consequence of (|3.40p that 

E \fn-{fn,gi,n)frdr<Gu;y2 I \fn-{fn,9i.n)frdr<GG' (3.50) 

SO imphes (IH^TZ)) ans thus 

Now, according to Proposition 13.31 and p.40p we have 

FUu^)>Xi,n'+2iTG f gt„,{r)rdr-GGLu-'/^ V |n-n*|2||/„||2 CG^. (3.51) 

Using and G < w^/^ we get 

-Fc.K)>Ai.„. +27rG / g^„.(r)rdr-GG2 
which was the missing lower bound to prove (|1.26p . Recall that G ^ uj^^'^ and g\ ^,(r)rdr oc o;-'^/^. 
S'iep 5. We prove pr27|) . We have 



Indeed 

Ai,„. +G/ |<?i,„.r >F^(u^) > Ai^„. +G / lu^r 

so 

|4 



/ \uX< I (3.53) 
Now, using the same calculations as those in the proof of Proposition 13.31 and (|3.40p we can prove 



>( gU,-G E ^"'/'l^-^*ni/nlli^(M+,.<i.)-C^G (3.54) 



Then plugging and ^^M^ into (153^ we get 



/ (KP-5L-)'<C E ^"'/'l«-^*ni/nlli2(K+.,rf,) + CG<G(G2c.-i/2 + G)<GG 
and the proof is complete. □ 
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4 The giant vortex state for a fixed coupling constant 



In this section we present the proofs of Theorems 11.21 and 11.31 They require improvements of the method 
presented in Sections 2 and 3. In particular we will need an energy expansion with a remainder term 
going to as w goes to infinity. We know from p.26|) that 

-FLK) = Ai,„. + 27rG / gt,,.rdr + 0(1) 

in the regime where G is a fixed constant. But, for any n, by the definition (|4.2p of 7„ 

Xi,n + 27rG / gt.n^dr > 7„ 

and one can realize from (|2.19p and (|4.7p below that the difference between these two quantities is of 
order 1. We thus need to analyze in some details the nonlinear problem defining 7„, and use the results 
to provide better upper and lower bounds to F^{ui^). We first present the analysis of the nonlinear one 
dimensional problems in Subsection 4.1 below. 



4.1 Nonlinear one dimensional problems 

We recall the definition of the one-dimensional energies 

EM) = 27r / (|/'(r)|2 + Vn{r)\f{r)\' + G\f{r)f) rdr (4.1) 

with 

jn = EMn)^inf\EM), f ^H\R+,rdr)nL'{R+,rVn{r)dr),27r f f{r)rdr^l}. (4.2) 



Note that is uniquely defined up to a multiplicative constant of modulus 1, that we fix by requiring 
that 

> 0. (4.3) 

We also define the rescaled functions p„ by 

p„(x) =c-i*„(i?„ + /i„a;) (4.4) 

where Cn is chosen so that ||/9„||i2 — 1: 

cl = - - . (4.5) 

27r/i„i?„ + 27r/i2 xpn{x)^dx 

We have 

Theorem 4.1 (Asymptotics for the ground states of the one-dimensional nonlinear problems). 
Let 7„ be defined by 

Suppose that n G A/1/2 o,^ defined in \2.11]) and let i?„ be defined in \2.4-^ . Let ^1 be the normalized ground 
state 0/ the harmonic oscillator \1.24-^ . 
We have as —> 00 and — > D 



7„ = K(i?„) + \l ^^^^ + I Uxfdx + J„ (4.6) 
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where J„ — 0{1). 

Let pn be defined by equation j^.^p . We have 

Pn--£,1+ hnTn + hlv^ + Oho(w"^/^) (4.7) 

where hn is defined in i2.12\) . Tn and ti„ are solutions to linear second order ODEs. They are bounded 
in harmonic oscillator norm, uniformly with respect to n. Also, for any p,q (z A/1/2 

Tp = Tg + OHo{^'^\p-q\) (4.8) 

From this theorem we deduce as a corollary 
Corollary 4.1 (Variations of 7„ with respect to n). 

Let n* and n* + I be defined as in Corollarv \2.1i Under the assumptions of Theorem \4-.l\ and for lo large 
enough, for any n n*, n* + 1 one has 

In = 7n- +C{n~ n*f (1 + 0(c^-i/2)^ (4 9) 

Moreover 

ln' + l>ln' -ClU-'^^. (4.10) 

We also state the equivalent of Corollary [ 



Corollary 4.2 (L^ bounds for Vl/n). 

For any n e A/1/2 a^'^ ^'^2/ 1 ^ P < +00 there is a constant Cp depending only on p such that 

ll*n|Up(R+,rdr) < CpL^'^^"^ (4.11) 

with the convention that ^ — if p ~ +00. 

The proofs are exactly similar to those of the corresponding results in Section 2, so we only give their 
main steps. 

Proof of Theorem \4-l\ We begin with an upper bound on 7„ using a test function of the form 

vTftest _ „test f I — Rr, 



hn 

We obtain 



7n<K(i?n) + \/-^^^^^+ , ^ „ I il{xYdx + 0{l) (4.12) 



2 2nhnRn 

We then prove an exponential decay property for and its first derivative exactly similar to that of 
Proposition l2.2l Step 1 of the proof of this Proposition needs no modification to apply to 
We use the maximum principle on the equation for which reads 

- A*„ + Vn{r)-^n + 2G*3 ^ (^^ + 2^^/+ |*n(r)|Vdr^ (4.13) 

Using (|4.12p we get easily 

-A*„ + {Vn{r) - Vn{Rn) ~ VK^R^) < 
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and we proceed as in the proof of Proposition 12.21 We then prove the equivalent of Lemma 12.31 which 
allows to get an expansion of 7„. We use this expansion as in Step 1 of the proof of Theorem 12. II and get 

for any e > 0, and follow the method of Step 2 of the same proof to improve the remainder term and get 

Pn=^l+OHo{0J-^^^+')- (4.14) 

We deduce an equation for p„ starting from (|4.13p : 

n" ^" n' I ^" 2 , O/^^a > 2 3 , l5 ^3 „ 

— Pn ~ J^P" 2 ^'-^c„n„p„ H - n^x Pn 



= hi (^7„ + 2nG J^^ \^nir)\^rdrj p„ + OL2{hl). (4.15) 

We then deduce from this equation that 

Pn = a + hnTn + /l,>„ + Oho{oJ'^^^)- (4.16) 

We have denoted t„ the solution of the problem 

and Vn the solution of the problem 

-u + X u — u — — jp/ina; ?i + 6 ""n-^ ^" 

where = 0(1) is chosen so that the right-hand side of the above equation is a function orthogonal to 
fi, which implies that the system indeed has a unique solution. In particular one has = J„ + 0(cj-i/2). 
This follows exactly Steps 3 and 4 of the proof of Theorem l2.1l Remark that the coefficients of the ODEs 
satisfied by t„ and Vn depend continuously on i?„ seen as a continuous variable for the range of n we 
consider. The esimates (j4.8p are a consequence of this fact. We deduce (|4.6p from (j4.7p . □ 

Proof of Corollary \4-l\ At this stage, we know that 

-fn=CiRn)+V{R„) + Jn 

where the function C is defined in equation (|2.75p . V is given by 

T^iRn) = 7T-l^ / ^i{^)'dx (4.17) 



27rhnRn 

and J„ = 0(1)- We have J„ — J„i 0{uj^^/^) for any m,n £ A/1/2, as a consequence of (|4.8p . and from 
the explicit expression of V one can compute that 

\V\R„)\<Clu'^' 

for any n E A/1/2 • We thus have, using Corollarv l2.1l 

7« = in'+C\n-n*\^ + 0{uj^/^)\Rn-Rn'\+Oiuj-^/^) 
= 7„- + C\n - n*f + 0{uj-^/^)\n -n*\ + 0{uj-^/^) 
= 7n-+C|n-n*|2(l + 0(c^-i/2)) 

for any n n* , n* + 1. 

Equation (j4.10p also follows from this discussion. □ 
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4.2 Some improved estimates 

In this subsection we aim at obtaining a better L°° bound for u^, which will be crucial in our analysis. 
Indeed, the first bound that we obtained (see Proposition l3.ip is far from being optimal as one can realize 
by taking the LP' norm of both sides of (jS.ip . The right-hand side has a norm proportional to w^''^ 
whereas ||uij||2^2(r2) — 1. A result of the analysis below is the improved bound 

|ka.||L»(R2) < Cuj^/\ (4.18) 

This will follow from (|4.20p and (|4.26p below and implies (this is a simple modification of the proof of 
Proposition 13. 1|) that (|3.ip can be improved as 

\uu{x)\ < Ccji/4e"'""ll="l"i|' for any x e E^. (4.19) 

The right-hand side of this inequality is now bounded in L^. 
We begin with a estimates for u^: 

Proposition 4.1 (Estimates for u^). 

There exist a constant a so that the following holds on 

\uuj{x)\ <Cwi/4g-aa.|b|-i|^_ (4 20) 

Moreover, for any < p < +oo 

with the convention that ^ — if p = -\-oo. 
Proof. We recall that 

uUr,0)-^ E (/n'5i,n)<?i,n(^)e"' 

so, using (|2.23p we get 

< J2 \{fn,gi,n)f\n-n*\A E 

< (4.22) 

where we have used (I3.46P to pass to the third line. The estimate (|4.20p is proved using the same kind 
of computations : We remark that if — 1| > Clj^^/^, then (|a;| — i?„)^ > (1 ~ C'C')(N| — 1)^ for some 
constant C. This is a consequence of the fact that if n € A/1/2 then i?„ — 1 = 0{uj^^^^). With a proper 
choice of C we thus have, using (j2.19p . that for any r S M so that |r — 1| > Cuj^^^^ and any n £ A/1/2 

\gi.n{r)\ < Cc^i/^e-^'*" I' < Cuj^/^e-"^^'^-^^". (4.23) 

We use this fact and the same trick as in (j4.22p to conclude that (|4.20p holds on the domain 1 1 |x| — 1 1 > Cto" 

On the complement of this domain e^'^'^''"^^' is bounded below, so (|4.20p is a consequence of (j4.22p . 
Finally fj^/U} follows from g^D])- □ 
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We now aim at improving p.40p . giving estimates in stronger norms. As already emphasized the 
most important result is (|4.26p . but we also state estimates (|4.24p and (|4.25p because they actually imply 
(|4.26p (see Step 5 of the proof below). 

Proposition 4.2 (Estimates in stronger norms). 

Recall that is defined in equation i3.9\) . The following estimates hold true when lo — > +00 and D: 

\\\uuj - u^IWh^r^) < CeOJ^ for any e > Q (4.24) 
\\\uu>-u^\\\m(m.^) < C^iv^^^ for any e > (4.25) 
\\uu, - u^\\l^{M2) < Ceuj" for any e > 0. (4.26) 

The proof uses mainly the Euler-Lagrange equation for m^, combined with elliptic estimates for the 
Ginzburg-Landau operator. We state the result we are going to use for convenience and refer to [38] or 
[3] for a proof. 

Lemma 4.1. Let A e (Vl^^'°°(M^))^ be a divergence-free map, g E L^(M^) and w a solution of 

- {V - iAf w = g inR^. (4.27) 
Then for any R > there exists a constant Cr > so that 
2 



V / \{d,-iA,)(dk-tAk)w\^ < 







([ 









(4.28) 



Moreover, Cr remains bounded as R goes to infinity. 
Proof of Proposition \J7^ Step 1. We claim that 

{u^ - u^) < C^^l^ (4.29) 

which implies 

lkc.-S^lki(E2) < (4.30) 

We begin by multiplying the Euler-Lagrange equation (|1.16p by u* . Integrating and injecting the result 
into {uuj — U(^), we obtain 



+ G {2\u^\^\u^\^ +u^\u^\^ul,+ul\u^\^u^-u^u*J -ulul) . (4.31) 
We note that a consequence of (j3.53p is 

\u^\^ <CGlo^/^. (4.32) 



The terms in the second line of (|4.3ip are also bounded by CGuj^^^ using Holder inequalities, (14.2ip and 
(|02)) . On the other hand, equations (fL26)) . ((06)l and (|4:20| yield 

F^uu) = F^{u^) + 0{Gu^/^) = Ai,„. + 0{Gu^/^). (4.33) 
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From (|3.5p and (|4.32p we also obtain 

= Ai,„. + 0(Gcc>i/2). (4.34) 

Then, using (|3.40p we deduce (|4.29p from (|4.3ip . remembering that in this section we assume that G is 
a fixed constant. 

Step 2. Using (j4.2f p together with (j4.30p and a Sobolev imbedding imphes 

\\u^\\Ln^-)<C,uj^/^+' (4.35) 

for any \ < p < +oo. We then interpolate between (remember that |l7i(^||j^2(][j2') = 1) and , make 
p +0O to get 

||w^||l«(R2) < for any e > 0. (4.36) 
Step 3. We now turn to the proof of (|4.25p . The Euler-Lagrange equation for gi ^ can be written 

2 

- (V - lujx^f g^^^e'^' + Dn'^ (jxp - l)%i,„e™^ = Al,„5l,ne"^ (4.37) 

SO that we get for 

2 

-{y-iux^fu^ + Dn'^{\x\''-lfu^= Ai,„(/„,5i,„)gi,„e*"^ (4.38) 
We substract this equation from the equation for u„ to obtain 

- (y - iuJX-^Y {u^ - U^) = fl^U^ - ^ \l,n if n, gi,n) gi,ne'''^^ 

,2 



Dn"^ (Ixp - 1)' (u^ - u^) - 2G\u^\''u^. (4.39) 



We denote 



S^:=H^u^- Ai,„(/„,5i,„)5i,„e"^-i:»OY(|xp-l)^(w^-u^)-2G|u„pw„ (4.40) 

and provide a bound in L^{B2r) to this quantity, for some R > 0. First, using (j3.40p and the exponential 
decay results of Propositions 13.11 and 14.11 we get 



Dlu;^ (|x|2 - 1)' - < C,cj'+^. (4.41) 



On the other hand 

2 



MwUw " ^ Ai,„ (/„,gi,„)5i,„e™^ 



C E (Al.n - Mt^)^ l(/n7.9l,ri) 
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We know that \^^\ < Cuj and — {tt^|1^2(R2) < Cuj -^/^ , whereas comparing (|1.26p and (|2.2ip we get 
|Ai,„ - < Cmax(w^/2, \n - Using ([5^ and recalhng that for n e 7Vi/2, |n - n*| < Cw^/^ 

we thus have 

2 

< Cw^lju^ - U,^||^2(K2) + ^ |(/«,5l,n)|^ 

neA/'i/2 

+ Cc^ ^ \n-n*\^\{fn,gi,n)\^ <Clj. (4.42) 
Gathering (|4.36p . (14.4ip and (|4.42p we have for any s > 

Applying Lemma 14.11 with A = ujx-^ wc find for any R 




(4.43) 



with Ce independent of R. Using the diamagnctic incquaUty (see [^) twice, we get 

\\V\u^~u^\\\m(BR)<C,u'+' (4.44) 
with Ce independent of R. This concludes the proof of (|4.25|) . 



E 

^6A^i/2 



(/n,5l,n) 9l,- 



Step 4- Let us now prove (|4.24p . We multiply (|4.39p by {u^ — Uui)* , integrate over and use the 
diamagnetic inequality to obtain 



< 2G 



(4.45) 

For both terms on the right-hand side of (I4.45P we use the Cauchy-Schwarz inequality and p.40p . combined 
with for the first one and for the second. This yields 

Step 5. We interpolate between and H^'^^ , make ^ and use a Sobolev imbedding to deduce 
(|4:26)) from (|4:24)) and g^S]). 

□ 



4.3 A refined lower bound on the interaction energy 

The last essential ingredient of our analysis is a refinement of Proposition 13.31 that we present in Propo- 
sition |43] below. Its proof uses the same ideas as that of Proposition 13.31 and requires to first improve 
the asymptotics for the modes /„, n G A/1/2 ■ We know from (|3.40p (G is now a fixed constant) that for 

n e Ni/2 

fn = (/„,ffi,„)ffi,„ + OL2(cc>-i/2). (4.46) 

We now have to be more precise and evaluate the term proportional to in the equation above. 

This is what we do in Proposition 14.31 The key observation is that (|1.27p and (|4.18p allow to linearize 
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the Euler-Lagrange equation for u^, at the price of a relatively small remainder term. Thus, inserting 
the Fourier expansion of wc obtain equations for each mode /„ that allow to improve (|4.46p . It then 
turns out that the term proportional to u!~^^'^ in (|4.46p has to be the solution of an elliptic problem that 
we describe in Lemma 

In this section we will often consider gi n and /„ (as well as any function defined on R+) as radial func- 
tions defined on K^. 

We denote by Ili.n and II^'-^ the othogonal projectors on the space spanned by gi,n and its orthogonal 
in H'^(R+,rdr) n L'^{R+, Vn{r)rdr), respectively. 

In the following Lemma we introduce some functions r„ that will be useful in the rest of the paper 
because they appear naturally when writing expansions for the modes (see Proposition 14.31 below). 
We also state the main properties of these functions that we will need in our analysis. In particular the 
estimate (|4.50p will be essential in the proof of Proposition [ 

Lemma 4.2 (Properties of r„). 
Let n G A/1/2- The problem 



'-Ar„ + VnTn - Ai,„r„ = n^„ {-2GLj^/^gl„,gi,n) 
/]g2 51,11^,1 = 0. 



(4.47) 



has a unique radial solution in H^{M.'^ ,rdr) C] L'^iM.'^ ,rVn{r)dr). It satisfies the bounds 

||r„|U2(R2) < C, ||r„||ffi(K2) < (4^48) 

and, for any e > 0, 

l|r„|U4(R2) <c,c^i/8+e^ ||r„|U6(R2) < ac^l/6+^ |ir„|Uoc(R2) < Cec^l/4+^ (4.49) 

Moreover, for any 71, m G A/1/2 o-i^d, any e > 

||r„ - r„|U2(R2) < CeW-i/2+=|n - m|. (4.50) 

Proof. It is easy to show that problem (14.47P has a unique radial solution in the energy space, for example 
by expanding both sides of the equation on the basis {gj.n)j=i...+oo- 
Using the expansion 

r« = (r™, 5j,n) gj,n (4.5i) 

combined with p.Sp the first line in (|4.47p can be written 

J2 ihn - Ai,„) (r„, 5,,n> 9i.n - i^-2GLu'/^ gl^. gi,„) ■ (4.52) 

Thus, taking the norm of both sides of (j4.52p we obtain 

YiXj,n - >^l,nf \{Tn,gj,n)f < II 5?,n- 5l.n 1 1 1,2 (r2) < Cuj'^ (4.53) 

because gf „. < Cw^/^ and gj^n is normalized in L^(M'^) for any j. Now, using (|2.18p and (|2.22p . one can 
find a constant < a < 1 such that 

Aj> - Ai,„ = oi^j.n + (1 - a)Aj,n - Ai,„ > a\j,n + (1 - a)A2,n - Ai,„ > aAj,„ > aAi,„ > Coj, (4.54) 
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for any j > 2. Thus, (|4.53p implies 

Y.\i^n,9j,n)f <C. (4.55) 

This yields the estimate in (|4.48l) because of the expansion (|4.5ip (recall that by definition (r„, = 
0). We also obtain from ((i35)) and 

thus the bound in (|4.49p holds. We have used again {Tn,gi.n) = 0. 
We now claim that r„ satisfies the pointwise estimates 



|r„(r)| <CeC^i/4+-e-*('^)' 
for some constant a, and for any a < 1 and r < Ca 

|r„(r)| <Ce-^-rS. 



(4.56) 



(4.57) 



The proofs are very similar to that of Proposition 12.21 so we give only the main ideas. First we remark 
that r„ is radial, so for any 77 > 

CUJ^/'^ > ||r„(a;)||//l(R2 = V2TT\\rn{r)\\Hi(R+,rdr) > \/27r7?|lr„(r)||^l([^_+oo[,<ir) 

and thus the bound 

follows by interpolation and Sobolev imbeddings as in Step 1 of the proof of Proposition [221 We deduce 
from (|iTf|) that |r„p satisfies 



- A|r„|2 + 2|r„|2 {Vn - Ai,„) < r„n^„ (-2Gl^i/X„.5i,») 

Using n^„(/) = / - ni,„(/) and „. < Cco^/^ we have 
n^„ (-2G^i/2^2^^^^^ 

Defining the domain 



(4.58) 



< C 



+ c 



1/2 2 2 



9i,n\ < Cugi, 



we thus deduce from (|4.58p that 



/„ = {xe M^ |r„| > 9i,n} 



-A|r„|2 + 2|r„|2(i4-Cc^) <o 



on /„. We prove that (|4.56|) and (|4.57p hold in /„ by using the maximum principle as in the proof of 
Proposition[121 On the other hand, and (H37)) hold on by definition, recaUing (PTTT)) and (H^H)- 

The LP bounds follow from 

We now turn to the proof of (|4.50p . We begin by noting that because of (|2.19p . a change of variables 
r — Rn + hnX and a Taylor expansion yield 



i||l2(r2) < Cuj ^^^\n - m\. 



(4.59) 



Thus 



R2 



^m9l,- 



Tm (ffl.m — 9l,n) 



< Cuj-^l'^\n-m\ 
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because ^mgi.m = and r„i is bounded in L^. This implies 



< Cuj~^/^\n-m\. 



On the other hand, denoting 
we get from (|4.47p the equation 

+ {Vm — Vn) T„i + {\l.m — ^l,n) 



(4.60) 
(4.61) 



+ iV„^ - Vn) r„. + (Ai,™ - Ai,„) r„,. (4.62) 



We estimate the norm of the right-hand side of this equation. First remark that, using the definition 
of the potential 14, we have for n,m & A/1/2 and |r — 1| < Cea;~^/^+^ 



\Vrn{r) ~ Vn{r)\ 



— 2uj(m — n) 



^\n~m\ m + n - 2u; + {m + n)0{uj 1/^+^) 

KGeOj^/^+'lm-nl. (4.63) 



We then write 



(K.-K)'r^ = 



|r~l|<Cw-i/2+e 



|r-l|>Cw-i/2+E 



(Kn-K)'r^ 



We estimate the first term using (j4.63p and the fact that Tm is bounded in L^. The second one is 
exponentially small because of (|4.56p and (|4.57p (recall that i?„ = 1 + 0{iL!^^^^) and /i„ cx w~^/^). The 
result is 

II {Vrn - K) r™|U2(R2) < Cet^l/2+e|„ _ (4 54) 

Also, using (j2.2ip . we have for n,m ^ A/1/2 



|Ai,™-Ai,„| <Cc^i/2|n-m|. 



(4.65) 



On the other hand 



L2(E2) 



2Gu;^/^gln,{gi,m-gi,n) <Guj^/^\n-m\ (4.66) 

because of (|4.59p . We now use the fact that n;j^„j — ^tn — ^i,n ~ ^i,m to obtain 

i^tm " ^i^i) (2Gw^^^ff?,„.ffi,n) = (^J^^ 2Ggln,gl „^ (gi,„ - ■^1'™''+ (^^^ ^^i'" ~ 9i,m)^ gi,r. 



and thus, using (|4.59p again, plus the Cauchy-Schwarz inequality for the second term 



Kn - nL) [2Guj^i-'gln'giA ^ ^^'^'1'^ - ™l- 



(4.67) 
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Arguing exactly as when we obtained (|4.53p from (|4.47p . the combination of (|4.62p and the estimates 
to (jiTBT)) yield 

J2 (^^■^" - Ai,„)' |(5„,„„ „)|' < CeLJ^+'\n - m\\ (4.68) 

Using (|4.54p we thus have 

$]|(5„,™,ff,,„)|' <C,a;-i+^|n-m|2 

and there only remains to recall (|4.60p to conclude the proof because {gj,n)j=i...+oa is a Hilbert basis for 
H^{M.+ ,rdr)rM'^{M.+ ,Vnrdr). ' □ 

We now show that the second order in (|4.46p is given, after an appropriate normalization, by the 
function r„ that we have introduced in Lemma 14.21 
Let us first introduce some notation: 

A/'i/4 = |n e Z, |n- w| < cji/^l (4.69) 

^74.+ = {^e AAi/4,(/„,.gi,„) ^0} (4.70) 
^72.+ = {^e AAi/2,(/„,gi,„) y^O} (4.71) 

In the rest of the paper we will always denote /?„ the n-th term of a generic sequence satisfying 

Y.Pl<C. (4.72) 

The actual value of the quantity /?„ may thus change from line to line. Note that /3„ can depend on oj 
and Da as long as (j4.72p is satisfied. 

Proposition 4.3 (Refined asymptotics for the mode /„). 

Let n e A/172, + • We have 

fn-{fn,gi,n)(9i.n+io-^^^Tn) < /3„ max (c^^V^, |n - n* ^c^-^/^A (4 73) 



and 



Fn [in ~ (/„,5l,n) (gi.n + ^'^'^V n)) < max (^"1/2, |n - 71* j^c^'^) (4.74) 

where r„ is defined in Lemma\ 



Proof of Proposition |^. j| We need an equation satisfied by /„ . First recall that (j4.20p and (|4.26p imply 
that 

lkc.|U»(R2) < Cuj^/\ (4.75) 
Then, using (|1.27p (G is now a fixed constant) we get 



This allows to rewrite the Euler-Lagrange equation for as follows: 

2 

- (V - iux^fu^ + {\^\^ - + '2Ggl„,u^ = fi^u^ + L, (4.76) 

where 

L = 2G {gl,^,u^ - \u^\^u^) (4.77) 
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satisfies H^Hl^ ^ Cu?-^'^. Using the Fourier expansion (|1.42p of we can write (|4.76p as 

^ e™^ (-A/„ + + 2Gff2 „./„) = ^ e"Vc./n + ^ (4.78) 

n n 

We multiply this equation by g-™^ and integrate over Q to get 

- A/„ + + 2Gg2 „./„ = ^^/„ + L„ (4.79) 

where 

L„(r)= / e-'^''L{r,d)dd, (4.80) 
Jo 

so that 

E ll^"lli^(R+,r.<ir) < < C^l/^. (4.81) 

n 

For n G A/1/2, + write /„ as 

/n = (/«, .gi,«) (<?!,„ + CJ-l/^Kn) (4.82) 

where 

(k«,5i,„)=0. (4.83) 

Remembering p.40p we have 

E l(/",5l,n)l'll'^n|li2(K+,,,,)<C. (4.84) 

We claim that 

II (/n, ffi.n) (k„ - r„) ||i2(R2) < /3„ max (uj'^l\ \n - n*\^^-'-') (4.85) 

and 

Fn ((/„,. gi,„) K - r„)) < /3„max (c^^/^, \n - n-fuj-"-^ (4.86) 

with /3„ satisfying (I4.72p . Let us denote 

C« := (/„,5i,n)K-r„) (4.87) 

Inserting (j4.82p into (|4.79p . using the fact that gi,„ is a normalized eigenfunction of — A + y„ for the 
eigenvalue Ai 

I 9l,n) l^n 

- 2Ggl^, {fn, gi,n) ^'^'^gi.n + (Mc. - Ai,„) ^j^^ ^^^^^ g^^^ ^ wi/2^„. (4.88) 
We multiply this equation by g\_n and integrate to get the useful relation 



(/«,5i.n) ( 2G / 5?,„.5?,„ - (/^o; - Ai,„) ) = / L„5i,„ - 2G(/„,5ri,„) / gi,„.K„gi,„ 

= 0(/?„c.3/4). (4.89) 

We have used Cauchy-Schwarz, the bound g\ n* < Gw^/^, (|4.8ip and (|4.84p . Substracting the equation 
for {fn,gi,n} Tn from (|4.88p . we obtain 

- ADn + VnDn - Ai,„Ai = (A^w - Ai,„) {fn,gi,n) ^ "^G {fn, gi,n) gj^n* l^n 

+ {fn, gi,n) ni,„ ([^J^^ - Ai,„) c.i/2<?i,„ - 2Ggl^,u^'^gi,,}j + c^^/^^n. (4.90) 
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Using (dUni), dsn and we obtain \fi^ - Ai,„| < Cmax(wi/2, - n*\'^) and thus 

IKm^^ - Ai,„) {fn,9i.n) Kn|li2(R2) < /3„ max (^w^/^, |?i - n*p^ (4.91) 
for any n G Ni/2.+- We have used (|4.84p and n* = lo + 0(1). Also, using 5? „. < Ct^^/^ and (gJH), 

pG {fn,gi^n) gin'^^n\\^,^^,) <^'/'/3„. (4.92) 

On the other hand (|4.89p imphes 

(/n,5l.n> ni,„ ((^io; " Ai,„) W^/^gi,„ - 2G'5i_„.LJ^/^5i,„ ) = 

w'/'(/n,5i,n) ((M..-Ai,„)-2G^^52„.g2^^ ll.9i.n|lL^(M2) < (4.93) 

We compute the norm of both sides of (|4.90p after having expanded the left-hand side on the basis 
(ffj".n)j"=i...+oo as in (|4.52p . Gathering equations (I4.9ip to (|4.93p . recalling (|4.8ip we obtain 

iKn - Knf \{Dn.9j,n)\' < C/J^ max (w^/^^ \n - n*|4) (4.94) 

where Aj^„ and gj^n are the eigenvalues and eigenfunctions of the restriction of — A + Vn to the space of 
radial functions. Recalling (|4.54p . (j4.94p implies 



which is (|4.85p because {Dn,gi.n) = (see the second equation in (|4.47p and (|4.83p V We also obtain 
from and 

F„(AO = hn\{Dn,gj,n)\^ < C ^ (Aj, „ - Ai,„) |(A.,5j.n)l' < Cplu;-^ max (^3/2^ |„ _ „*|4^ ^ 

J>2 i>2 

thus (|4.86p holds. We have used again {Dn,gi,n) =0. □ 

We are now able to present the main result of this section, a refinement of the lower bound on the 
interaction energy. Its proof uses the same as ideas as that of Proposition 13.31 but the next order in the 
lower bound is now computed. This is made possible because of the refined asymptotics of Proposition 



Note however that the right-hand side of (|4.73p increases with \n — n*\ so the refinement is not equally 
efficient for all modes. On the other hand p.46p states that the larger |n — n*| is, the smaller is the 
mass of the corresponding mode, so the smaller is its contribution to the interaction energy. It turns 
out that the best lower bound available with our results can be computed by using the refined results of 
Proposition 14. 31 onlv when n G A/1/4.+. 

Proposition 4.4 (Refined lower bound on the interaction energy). 
We have, for any e > 0, 

/ \uX> E \ifn:9l.n)H (gi.n + UJ-'^^tS rdr 

+ E \{fn:gi.n)f f g^vdr ~ C .^'^ ' . (4.95) 
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Proof. We begin by using Lemma l3.1l 

/ \uX>2n I \fp\^\U\'^dr (4.96) 

and recall that we can restrict the integration domains to 1/2 < r < 3/2 at the price of an exponentially 
small remainder (see p.20p ). We implicitly do this restriction in the rest of the proof. 
We continue to denote /3„ the n-th term of a generic sequence satisfying (|4.72p and note that for any 
n G A/i/2 we have 

ll/n||L2([l/2,3/2],dr) < C"! I /n || (K+ .rdr) < Pn 

because of the LP' normalization of . Moreover 

ll/n||Hi([l/2,3/2],dr) < C| I /n || ffi (R+ ,rdr) < Pnl^^^^ 

follows from the bounds J2n ^n{fn) < Fuj{uuj) < Civ. Interpolating between these two estimates we get 

ll/n|lHi/2([i/2,3/2],dr) < Pn^^^^^- 

Then we obtain, for any e > 

ll/nlUa([i/2,3/2],d.) < C,/3„c^i/s+^ (4.97) 

by the Sobolev imbedding of H^^^ in for any p and interpolation between and with p +oo. 
If n ^ A/i/4 this estimate is improved: 

ll/n|U«([l/2.3/2],dr) < C,l3nL^~'/'^+' . (4.98) 

Indeed, for n ^ A/1/4 we have 

ll/n|lL2([l/2,3/2],dr) < C\\ fn\\ (R+ ,rdr) < PnUJ~^^'^ (4.99) 

using (|3.46p . Moreover the expansion 

i^«(/n) = E^J-."l(/"'ffj-.")l' (4.100) 

combined with (piiS)) . |(/„,5i^„)|^ < /3,V"^^^ and Ai,„ = 0(cj) yields 

ll/n||Hi([l/2,3/2],<ir) < C\\ fn\\ m (R+ ^rdr) < C \/FjJn) < Pn^J^^^- (4.101) 

Interpolating between (|4.99p and (j4.10ip we obtain (|4.98p . We note that there also holds, for any n E A/1/2, 

II (/n,5l,n)5l,n||Le([i/2,3/2],dr) < C/^nw'^' (4.102) 

and for any n </ A/1/4, + 

II (/«,<?l,n)5l,n|lL«([l/2,3/2],<ir) < CPnUJ-^/^^. (4.103) 

These last two estimates are basic consequences of the properties of gi^n, see Corollary[ 
We now bound from below the p-th term of the sum (j4.96p . distinguishing two cases: 
Case 1: p ^ 

We have, for any g e A/1/2 



\fp?\U? - l(/p>5i,P>h(/9:5i,9>l^ffi,p5i,g) ^'^^ 



< 



C'Wfq - {fq,9l.q) 9Uq\\L4fp\\h (II/«IIl6 + II {fq.9l.q) ffl^glUe) 

+ C\\fp- (/p,5l,p)5l,plUHI (/g>5l,g)ffl,?lli6 (ll/plU« + II (/p,5l,p>5l,glU«) 

< Cio-^'^+'plPl (4.104) 
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using p.40p and the estimates (|4.97p to (|4.103p . Wc then use exactly the same technique as in the proof 
of Proposition 13 . 31 to obtain that for any q G A/1/2 

/ 9lr>9l,rdr > J gl^rdr - Cu;-'/'\p - - C. (4.105) 

We conclude from (|4.104p and (|4.105p that for any p ^aAi/4.+ 

E ^^J\fpWrdr> J2 \{fp,9i,p)f\{f,,9i,,)f(^27:Jgl^rdr-Cu;-'/'\p-q\'-c'^ 

- C,l3luj~^'^+'. (4.106) 



Case 2: p e A/i/4_-|_ We again distinguish two cases. 

Case 2.1: q ^ Ni/i,+ Wc argue as in Case 1. Using estimates (|4.97l) to (|4.103p . 

{\fp?\fr,? - \{fp,9i,p)?\{U9i.,)\'9lp9l,) rdr < C,u;-'/'+' (3^^ 
because q ^ A/i/4_+. Also 

/ 9lp9l,rdr > J gl^rdr - Clu~'/'\p - q\^ - C. 
and it is clear from (|2.19p and the estimates in Lemma that for any e > 

glprdr^ J (51^^ + c^-i/%)%dr + O(w^). 
Gathering (|4.107p . (|4.108p and (|4.109p we thus have for any p S Afi/4,+ 

J2 2n I \fp\^\f,\\dr > \{fp,9i,p)\'\{U9i,,)f I [g^^p + cv-'/^TpYrdr 

Case 2.2: q G J^i/4,+ This is the case where the refined asymptotics of Proposition will be crucial. 
Note that for n £ A/1/4,4. ^4:.73^ becomes 



(4.107) 



(4.108) 



(4.109) 



(4.110) 



We begin with 

\fpW - \{fp,9i..p)f\{U9i,,)\' [gLp + Lu-'^'TpY (.91,, 



(4.111) 



-l/2p 



c 



rdr < 

h-{h,9i..,) (5i,,+^"'/'r,)||^J|/p||ie(||/,|Ue + ||(/,,5i,,) (gi,, + c.-i/2r,)||^J 



c 



fp - {fp,9i,p) (si.p + ' 



{fq,9l,q) 



-l/2p 



L6 



(||/pIIl« + \\{fp,9i,p) + < C,u;-'/'+'(3X- (4-112) 
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We have used (|4.97p and (|4.102p along with (|4.11ip and the estimates of Lemma li?^ 

Now we use the bounds llgi.pIlL* < Cuj'^/^ (see Corollary 12. 2p and ||rp||i4 < CeCJ^/^^^ (see Lemma[ 

to obtain 



.91. P + ^ 



-l/2p 



+ 2lu-'/^ J gl^Tpgi^prdr - C,u;~'/^+'. (4.113) 



Using (|2.14p and (|2.19p we compute 



J gi.pgi.qrdr > hpCj^pC^glp.g 
with a change of variables r = Rp + hpX where 

Ip.q = Jp,q + "^Kp^q + 2Lp^g + 0{uj-^) 

and 



J, 



p,q 



Rp Rq 



Rv — R, 



{Rp + hpx) dx 



^p,q 



h„ hn 
hp 
h„ 



{Rp + hpx) dx. 



We recall that for any p,q A/1/2 
and 



\hp - hq\ fxui - q\ 



\Rp ~ Rq\ ecu ^\p~ q\. 
Then arguing exactly as in the proof of Proposition 13 . 31 we get 

Jp,q = Jp.p + 0{uj^'^\p - gp) 

whereas simple Taylor inequalities coupled with the estimates (|2.20p yield 

Kp,q = Kp^p + 0{uj^^\p - q\) Lp^q = Lp^p + 0{uj^'^\p ~ q\). 

Moreover, (|2.19p and (|2.20p imply that for any p,q e 

\clp-clj<Cu;-'\p-q\. 

All in all, 

2n I glpglqvdr >2nJ g^p - Cu'^'\ - qf - Cu;-'/'\p - q\ - Cuj' 



1/2 



On the other hand 



J 9lp^q9l,q- J glp^P < {jst-X)^ {J 



1/2 

{ghq - gi..pY ) + ( / gti 



1/2 



(4.114) 



(4.115) 



) {Rp + hpx) dx (4.116) 



(4.117) 



(4.118) 



1/2 



< \\9i.p\\\<i\\^q\\LA\9i.p - 9i,q\\L^ + llfi.pllielirp - rq||i2 

< C,u'\p~q\ 



(4.119) 
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using Corollary [221 Km . g^Ol) and (11351) . 

We conclude from (|11T^ . (jilT^ . (jilTHl) and that for any p e A/'i/4,+ 



4 

rdr 



E 2n f \fp\^\f,\\dr > \{fp,gi.,p)\'\{U9i..,)f (in [ [gi,p + u;-'/% 

-36AAi/4,+ 9eA/'i/4,+ ^ 

To conclude the proof we sum (j4.106p , (|4.110p and (|4.120p and use the following estimates that we 
obtain from ((OOl) and ((Oe)) : 



E l(/p-5i,p)l'l(/„5M)l'b-# < E l(/p,5i,p>l'l(/,,5i,,>l'b- 



* |2 



+ E \(fp'9hp)f\{fq,gi.q)f\q 

P-qeN'l/2 

and 

E \Up,9l,p)?\{fq.9l,q)?\p-q\ < I E \\fp\\\HR+,rdr)\\fqfL^(K,+ ,rdr)\p- 

\ 1/2 

E \\fp\\L'^(R+,rdr)\\fq\\'L^(VL+,rdr) j ^ C*. 
\p,q&Ml/2 I 



X 



Also, 

E \{fq.9l,q)?<Cu^''''- 

g^A/'i/4,+ 

Finally 



□ 



4.4 Proofs of Theorems Q and Q 

Proof of U.29]) . The upper bound is obtained by taking V&n.e™ ^ as a test function in F^^. For the lower 
bound we first write, using Proposition 14.41 

> E^"(/")+ E 2^G|(/„,.gi,„)|' / (<?!,„ +cj-i/'r„)%dr 

E 27rG|(/„,gi,„)|' / 5t„rdr - CeC.-l/4+^ (4.121) 



"6A^i/2\A^i/4.- 
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Now, using Lemma [^?T1 (|4.74p and the definitions of Ai.„, A2,ti and En we obtain from (|4.12ip 

,n + W '/'r„) rdr 



+ 



l(/n,5i,n)|' (Fnigi,n) + 2nG [ gl„rd,) + [ 



R+ 



> Y l(/-5l.n)|'^n(5l.„+C^-'/'r„ 



ll/n|lL2(R+,r£ir)^n (5l,n) 

neA^l/2\M/4,+ 



(4.122) 



where c is the constant appearing in Lemma l2.1l We have used the following : 

A2,„ > Ai,„ + Cu;> En{gi,n) + Cu - C'lo^'^ > S„(gi,„) 

which is proved by computing En{gi.n) with the results of Theorem 12.11 and Proposition 12.11 in mind. 
The normalization of gi^n and the orthogonality of gi^„ and r„ imply 



27r 



L^{R+.rdr) 



> 1 



SO 



En [gi,n+iO-^^''Tn) > 2tT 

On the other hand, because of (|4.73p we have 

27i'll/«|li2(R+_rd^) = 27r |(/„,5i,„)|^ 



gi.n + 



gi,n + UJ' 



2,-3/2 



■/3> 



(4.123) 



(4.124) 



for any n £ Afiu,+. Using that 7„ = 0(w) we finally obtain from (I4.122p . (I4.123P and (|4.124p 



''1/4,+ 

-FL;(Uw) > E 2'''ll/nllL2(R+,r£ir 



)7„+ E 27r||/„| 



L2(R+^rdr) 



(7„. + Cc^) - C,u-^'^+'. (4.125) 



Note that for n ^ A/1/2 one has Ai_„ > cw > 7„. + Cuj because c > \/6 and 7„.« < -v/Gw which can be 
computed from (|4.6p . The lower bound in (|1.29p follows from (j4.125p and Corollarv l4.1l recalling 



□ 



Proof of m^B- Using 



we obtain from (|4.125p 



rdr) 



ISA/", 



1/2 



1/2 
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with C > 0, thus, Corollary 14. II implies 

C,^-V4+.> ^ Mfn\\l^M+,rdr)\r^-rir+ E Mfn\\hiM+..rdr)C^ 

This immediately yields 

E 27r|l/„|li.(«+_,,,)|n-n*p < ac.-i/4+^ (4.126) 

and 

E 27r||/„||i.(«,_,,,)<C,c.-^/4+^ (4.127) 

We now introduce 

= a«'„.e"*« + /3*„-+ie'("*+i)^ (4.128) 

where 



a := v27r||/„. ||i2(R+ rdr) 

/? := V2^||/n- + l||L^(K+,rdr)- (4.129) 

We note that and (gT]) imply 

gi.n = + Ol2(r2) (w-i/2) (4.130) 
so, using equation (j4.126p and (|4.127p combined with (|4.73p we see that after extraction of a subsequence 

\\u^ - CUa,||L2(R2) < CeUJ^^^^^^ 

with c a constant of modulus 1, so that we have, recalling (|4.1ip and (|4.36p 



\Uuj\'^ - \Uuj 



R2 



I"*! < C" (||Mc^||i6(R2) + ||Ut^||i6(R2)) \\U^ - CUu;\\l^(R2) 

< - «i^|U2(R2) < CeC^3/'+^ (4.131) 

The estimate (|4.7p yields by a change of variables similar to those we used in the proofs of Propositions 
landOl 

/ \^nA^=l |*„. + l|^ + 0(c^-l/2)- / |vl/„.|2|vl;„.^^|2 + 0(^-1/2) 

JR2 Jr2 Jr2 

Thus, computing the interaction energy of ii^, using (I4.13ip and the fact that + (3"^ — \ + 0(aj~^/^+^) 
for any e > we have 

In' > F^[u^) > 7„. + 2a2/32 f |^„. |4 _ ^^1^^/**+^ 

JR2 

JR2 

But, according to our results (|4.126p and (|4.127p . either or f3'^ is bounded below. Then, along some 
subsequence we have either < CeW~^/^+^ or < CgO;"^/®"'"^. Renaming n* if necessary, we have 
proved that along some subsequence 

E \\fn\\lHM^,rdr) < C^^-'^'^' (4.132) 

which concludes the proof of (|1.30p . recalling (|4.73p and (|4.130p . 

□ 
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Finally, we present the 
Proof of Theorem \1.3[ Let us define 

A/'7/i6,+ := {« e Z, |n -n* I < 0.7/16, (/„,gi,„) ^O}. (4.133) 

and 

u^r.e):^ J2 {fn,gi,n)(bne"'' (4.134) 

"£.^7/16.+ 

where we have denoted 

0„ :=5i,„+c^-i/2r„ (4.135) 

for short. We claim that 
Firstly we obtain 



L2(R2) < ^ 27r||/„||^2 + ^ Wfn- {fn,gi,n/ Vn\\L2 



from (|4.126p and (|4.127p . We then note that satisfies 

- A(f>n + Vn<j>n = Ai,„0„ + H^,, (-2Gw^ gl,,, gi^n) ■ (4.138) 
Multiplying (|4.138p by {fn,gi,n) e*"^ and summing over n G J^7/ie,+ we deduce that 



- {V - iLux-^y + V{r)u^ = Ai,„ (/„,gi,„) 0„e"* 

+ E (/n,5i,™)n^„(-2Gc^i/Xn*5i,n)e"^ (4.139) 

ne7V7/i6,+ 

Substracting this equation from the Euler-Lagrange equation p.l6p we get an equation for — u^: 

- (V - iujx-^)" {u^ - u^) + V{r) {u^ ~ u^) = /i^Uc^ - 2G\u^\'^u^ 

- E Ai,„(/„,5i,„)0„e"^- J2 (/n,5i,n>nt„(-2G^i/V,„-5i,n)e'"^ (4.140) 

ne7V7/i6,+ neA/'7/i6,+ 

The same technique as in Steps 3, 4 and 5 of the proof of Proposition 14.21 allows to deduce (|4.136p from 
(Q-TTD and (Oim . 

Proof of Item 1 . 

We write, using (|4.136p 

||u^|-|<^n.||<|l-|(/n*,5l.n.)|||<^„.|+ E , 5l,„) | | + C,C.-^/''+% (4.141) 

A consequence of Lemma W?]\ (more precisely the L°° estimate in (|4.49l) ). (|2.19p and (|4.7p is 

<^„ = + Ol- (c."^/*) = *„ + Oloc (c^-1/4) = c„e-'*('^)' + Oioo (w-V4) (4.142) 
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where c„ is defined by (|4.4p and d is the standard deviation of fi. A consequence of (14. 7p is 

because is an even function. We recall that, for some constant h 

Rn = l + o(cj-i/2) (4.144) 

h,, = hLj-^/^{l + o{l)) (4.145) 

for any n E J^j/ie, so arguing as in the proof of Proposition 14. II one deduce from (|4.14ip that 

\\U^\- <cl\l- !(/„., II + l(/-5l,n)l J 

+ CeW"^/^2+^ (4.146) 

with 

. := (4.147) 

On the other hand (|4.132p and imply 

\l-\{fn',9l,n')\\<C,iU-^/^+' 

and 

E l(/".5l,n)l < C- E \\fn\\hiM^,r,r)\n-n*\A E l^-^l"' 

so we obtain (fT3T1) from (|4.146p . 

Proof of Rents 2 and 3. These final results are consequences of (|1.3ip . Using the facts that 

= c„.e~''(^)' +Oloo(c^-i/4) 
and Cn* oc w^^^ one can see from (jl.3ip that Ui^{xuj) uniformly on any region where 

II^E^I - 1| > Scj-Hn{u;) 
with S > a. On the other hand, if = r^^e^^"^ is a point such that 

Uuj{x^) 

as w — > oo, (|1.3ip impfies 

for any constant C as r r^^ and thus 

K - -Rn- 1' > ^K' (In(^) - C) > a^-i(ln(L^) - C) (4.148) 
which completes the proof because 

Rn' ^1 + 0{lO-^). 

□ 
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